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Abstract 

We have studied hadronic events from e+e“ annihilation data at centre-of-mass energies 
from 91 to 209 GeV. We present distributions of event shape observables and their mo¬ 
ments at each energy and compare with QCD Monte Carlo models. From the event shape 
distributions we extract the strong coupling Og and test its evolution with energy scale. 
The results are consistent with the running of Og expected from QCD. Combining all 
data, the value of ag(Mz) is determined to be 

ag(Mz) = 0.1191 ± 0.0005 (stat.) ± 0.0010 (expt.) ± 0.0011 (hadr.) ± 0.0044 (theo.) . 

The energy evolution of the moments is also used to determine a value of Ug with slightly 
larger errors: as{Mz) = 0.1223±0.0005(stat.)±0.0014(expt.)±0.0016(hadr.)^Q;QQ3g(theo.). 
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1 Introduction 


Hadronic final states produced in the process e’''e“ —qq are a valuable testing ground for 
the theory of the strong interaction in the Standard Model, Quantum Chromodynamics 
(QCD). The hadronic system in the energy range considered here is complex, consisting 
of typically 20-50 hadrons. Many “event shape” observables have been devised which 
provide a convenient way of characterizing the main features of such events. Analytic 
QCD predictions of the distributions of several of these event shape observables have 
been presented in the literature (see e.g. ref. [1]), and can be used to determine the cru¬ 
cial free parameter of QCD — the coupling strength ttg. These predictions describe the 
distributions of quarks and gluons, while the distributions of hadrons are measured in 
the data. In confronting the data with theory, Monte Carlo models of the hadronization 
process are commonly used to relate the partons and hadrons. Analytic QCD predictions 
have also been made for the moments of event shape distributions, whose evolution with 
centre-of-mass (c.m.) energy permit complementary determinations of ag. The determi¬ 
nation of tts from many different observables provides an important test of the consistency 
of QCD. In addition, measurements of event shape distributions have proved invaluable 
for testing and tuning Monte Carlo models of hadron production in e+e“ —hadrons. 

In this paper we present a coherent analysis of event shape distributions and moments 
using data collected by the OPAL detector at 12 c.m. energy points covering the LEP c.m. 
energy range of ^/s = Ec.m. = 91-209 GeV. Results at 192-209 GeV are published for the 
Erst time here. Partial results at 91-189 GeV have been published by OPAL previously [2- 
5]; these are superseded by the present measurements in order that the data at all energies 
can be analysed and interpreted in a consistent manner. The results presented here use 
identical analysis procedures throughout, and some event shape observables are included 
for the first time. In several cases, improved theoretical calculations are now available, as 
described in Sect. 4. The results at 91 GeV are based on calibration data taken during 
the LEP II running period (the period from 1996 onwards when LEP operated well above 
the Z mass); these share the same detector configuration (slightly different from that used 
in the earlier LEP I phase when LEP operated close to the Z peak) and reconstruction 
code as the higher energy data, which means that we can compare results over a wide 
energy range with minimal systematic differences between energies. Similar results from 
other LEP collaborations can be found in refs. [6-15]. Another recent OPAL paper [16] 
uses the same data sample as the present study to measure jet rates. 

The structure of the paper is as follows. In Sect. 2 we give a brief description of the 
OPAL detector, and in Sect. 3 we summarize the data and Monte Garlo samples used. 
The theoretical background to the work is outlined in Sect. 4. The experimental analysis 
techniques are explained in Sect. 5 before the measurements are presented and compared 
with theory in Sect. 6. 


2 The OPAL detector 

The OPAL detector was operated at the LEP e’''e“ collider at GERN from 1989 to 2000. 
A detailed description can be found in ref. [17]. The analysis presented here relies mainly 
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on the measurements of momenta and directions of charged particles in the tracking 
chambers and of energy deposited in the electromagnetic calorimeters of the detector. 

All tracking systems were located inside a solenoidal magnet which provided a uniform 
axial magnetic held of 0.435 T along the beam axis^. The magnet was surrounded by a lead 
glass electromagnetic calorimeter and a hadron calorimeter of the sampling type. Outside 
the hadron calorimeter, the detector was surrounded by a system of muon chambers. 
There were similar layers of detectors in the forward and backward endcaps. 

The main tracking detector was the central jet chamber. This device was approxi¬ 
mately 4 m long and had an outer radius of about 1.85 m. It had 24 sectors with radial 
planes of 159 sense wires spaced by 1 cm. The electromagnetic calorimeters in the barrel 
and the endcap sections of the detector consisted of 11704 lead glass blocks with a depth 
of 24.6 radiation lengths in the barrel and more than 22 radiation lengths in the endcaps. 


3 Data and Monte Carlo samples 

The data used here were recorded from 1995-2000 using the OPAL detector at LEP. 
In 1995 the LEP c.m. energy was increased above the vicinity of the Z peak in runs at 
T'c.m. = 130 and 136 GeV. By 2000, the maximum c.m. energy had reached 209 GeV. All 
of the data recorded above the Z peak are analysed in the present study. In addition, 
interspersed at various points during the high energy LEP running, calibration runs were 
taken on the Z peak, at ^/s = 91.3 GeV. These data were recorded with identical detector 
conhguration and performance, and reconstructed with the same code, as the high energy 
data. For the purpose of analysis, the data have been grouped into small energy ranges, 
which are often merged into larger ranges for clarity of presentation. Table 1 summarizes 
the c.m. energy points used, the integrated luminosities at each point and the numbers of 
events employed for analysis after the selection described in Sect. 5.1. 

Samples of Monte Garlo simulated events were used to correct the data for experi¬ 
mental acceptance, efficiency and backgrounds. The process e’''e“ —qq was simulated 
using JETSET 7.4 [18] at ^/s = 91.2 GeV, and at higher energies using /C/C2f 4.01 or 
/C/C2f 4.13 [19] with fragmentation performed using PYTHIA 6.150 or PYTHIA 6.158 [18]. 
Gorresponding samples using HERWIG 6.2 [20] or /C/C2f with HERWIG 6.2 fragmenta¬ 
tion were used for systematic checks. Four-fermion background processes were simulated 
using grc4f 2.1 [21] or KORALW 1.42 [22] with grc4f [21] matrix elements and with frag¬ 
mentation performed using PYTHIA. The above samples, generated at each energy point 
studied, were processed through a full simulation of the OPAL detector [23], and recon¬ 
structed in the same way as real data. In addition, for comparisons with the corrected 
data, and when correcting for the effects of fragmentation, large samples of generator- 
level Monte Garlo events were employed, using the parton shower models PYTHIA 6.158, 
HERWIG 6.2 and ARIADNE 4.11 [24]. Each of these fragmentation models contains 
a number of tunable parameters; these were adjusted by tuning to previously published 

^In the OPAL coordinate system the x-axis points towards the centre of the LEP ring, the t/-axis 
points approximately upwards and the z-axis points in the direction of the electron beam. The polar 
angle 9 and the azimuthal angle (j) are defined w.r.t. z and x, respectively, while r is the distance from 
the z-axis. 
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OPAL data at y/s ~ 91 GeV as described in ref. [25] for PYTHIA/JETSET and in ref. [26] 
for HERWIG and ARIADNE. 


4 Theoretical background 


4.1 Event shape observables 

The properties of hadronic events may be described by a set of event shape observables. 
These may be used to characterize the distribution of particles in an event as “pencil-like”, 
planar, spherical, etc. They can be computed either using the measured charged particles 
and calorimeter clusters, or using the true hadrons or partons in simulated events. The 
following event shapes are considered here: 


Thrust T: dehned by the expression [27,28] 


T 


max 

n 


fTn\Pi:^\ 

V Y.^\p^\ ) 


( 1 ) 


where pi is the three-momentum of particle i and the summation runs over all 
particles, which may be the measured particles, or the true hadrons or partons in 
Monte Garlo events. The thrust axis ut is the direction n which maximises the 
expression in parentheses. A plane through the origin and perpendicular to 
divides the event into two hemispheres Hi and H 2 . 


Thrust major T^aj.: The maximization in equation (1) is performed subject to the 
constraint that n must lie in the plane perpendicular to ut- The resulting vector is 
called Ut ■ ■ 

Thrust minor Tmin.: The expression in parentheses in equation (1) is evaluated for the 
vector Ut ■ which is perpendicular to both Ht and Ht ■ ■ 

Oblateness O: This observable is dehned by O = Tmaj. — Tm^n [29]. 

Sphericity S and Aplanarity A: These observables are based on the momentum ten¬ 
sor 

= «,/?= 1,2,3, (2) 

zZiPi 

where the sum runs over particles, i, and a and (3 denote the cartesian coordinates 
of the momentum vector. The three eigenvalues Qj of S°‘^ are ordered such that 
Qi < Q 2 < Qs- These then dehne S [30,31] and A [32] by 

S = ^{Qi + Q2) and A=^Qi. (3) 


C- and O-parameters: The momentum tensor is modihed to become 


oa/3 _ J2iiP?Pi)/\Pi 

Ei \Pi\ 


a, /5 = 1, 2, 3 . 


( 4 ) 
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The three eigenvalues Xj of this tensor dehne C [33] through 


C — 3 (AiA 2 + A2A3 + A3A1) 


( 5 ) 


and D through 


D = 27 A 1 A 2 A 3 . 


( 6 ) 


Jet Masses Mu and M^,: The hemisphere invariant masses are calculated using the 
particles in the two hemispheres Hi and H 2 . We dehne Mh [34,35] as the heavier 
mass, divided by i/s, and Ml as the lighter mass, likewise divided by i/s. 


Jet Broadening observables Bx? 

quantity 


Bk 


Bn and Bw: These are dehned by computing the 


\Pi X ^t|A 

V 2E.b*l J 


( 7 ) 


for each of the two event hemispheres, Hk, dehned above. The three observables [36] 
are dehned by 


Bt = Bi + B 2 , Bn = min(Bi, B 2 ) and Bw = max(Bi, B 2 ) ( 8 ) 


where Bx is the total. Bn is the narrow and Bw is the wide jet broadening. 

Transition value between 2 and 3 jets The value of the jet resolution parame¬ 
ter, I/cut, at which the event makes a transition between a 2-jet and a 3-jet assign¬ 
ment, for the Durham jet hnding scheme [37]. 

In the following discussion, whenever we wish to refer to a generic event shape observ¬ 
able we use the symbol y. In almost all cases, larger values of y indicate regions dominated 
by the radiation of hard gluons and small values of y indicate the region inhuenced by 
multiple soft gluon radiation. Note that thrust T forms an exception to this rule, as the 
value of T reaches unity for events consisting of two collimated back-to-back jets. We 
therefore use y = 1 — T instead. 

For all of these event shapes, a perfectly collimated (“pencil-like”) two-jet hnal state 
will have y = 0. 0{as) QCD processes generate planar qqg conhgurations; for most of the 
observables, these processes will generate contributions at y 7 ^ 0 — these are sometimes 
referred to as “three-jet” observables. However, hve of the observables (Tmin,, Ml, Bn, D 
and A), are still zero at 0{as), i.e. for planar events, and receive their leading contributions 
at 0(a'^) — these are referred to as “four-jet” observables. 


4.2 QCD predictions for event shape distributions 

QCD perturbation theory may be used to make predictions for event shape observables [1]. 
In order that these predictions be reliable, it is necessary that the value of the observable 
be infra-red stable (i.e. unaltered under the emission of soft gluons) and collinear stable 
(i.e. unaltered under collinear parton branchings). 

The QCD matrix elements in e’''e“ annihilations are fully known to 0{ag) [38], i.e. 
to next-to-leading order (NLO) for those observables dominated by three parton hnal 



states. In the two-jet (low y) region, however, the effect of soft and collinear emissions 
introduces large logarithmic effects depending on L = log(l/|/), such that the leading 
dependence on ag and L at each order, n, is proportional to For six of the three- 

jet observables studied here ((1 — T), Mh, i?T, .Bw, C and these large logarithms 
can be resummed to next-to-leading order, referred to as the next-to-leading-logarithmic 
approximation (NLLA). The most complete QCD predictions come from combining the 
0 ( 0 ;^) and NLLA predictions, taking care not to double count those terms which are in 
common between them. Further details may be found in ref. [2,39], and only a brief 
outline of the procedure is given below. 

The QCD calculations make predictions for the cumulative cross-section 

r 1 da 

Jo 

which take the following form for the NLLA calculations [40-43]: 

RNLLAiy) = (1 + C{ag -h Csttg) e-x^[Lgi{agL) -h g2{asL)] , (10) 


where for brevity we write ttg for {ag/2T[). The corresponding formula for the O(a^) 
prediction is: 

Ro{ai){y) = ^ + A{y)ag + B{y)al . (11) 

For the analysis presented here, the “log(i?)” matching scheme [40] is adopted for 
combining the Ola"^) and NLLA predictions. This matching scheme involves taking the 
logarithm of equation (11) and expanding as a power series, yielding: 


^^Ro{al){y) = A{y)ag + 


m 



al + 0{al) 


( 12 ) 


and similarly rewriting equation (10) as: 

lni?NLLA(|/) = Lgi(asL) g2{oigL) -\- Ciag [(F 2 — —C^]ttg -|- 0{ag) . (13) 

In the log(i?) matching scheme the terms up to 0{a‘^) in the NLLA expression, equa¬ 
tion (13), are replaced by the O^a^) terms from equation (12). 

Both the 0{ag) and NLLA QCD predictions depend on the choice of renormaliza¬ 
tion scale, fi (see ref. [1] for example). The renormalization scale factor is dehned as 
Xfj, = y/Ec.ra. where as{y) is the expansion parameter which appears in the NLO per¬ 
turbative predictions above. Naively y would be expected to be of order of, but not 
necessarily equal to, A QCD calculation to all orders should be independent of 

but a truncated fixed order calculation does in general depend on For exam¬ 
ple, in the 0(ag) calculation, the second order coefficient jl3(y) has to be replaced by 
jt3(y) + Pq Inx^Aiy) where /do = H — \n-p is the leading order /d-function coefficient of the 
renormalization group equation and np = 5 is the number of active quark flavours. Similar 
modihcations apply to the NLLA calculations. In addition, the theoretical cross-sections 
are usually normalized to the Born cross-section ao while the experimental distributions 
are normalized to the total hadronic cross-section, atot., which itself depends on ag- 

ffiot. = o'o (^1 + (14) 
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and this is taken into account by means of the replacement Biy) B{y) — 2A{y). 

The analysis of cts presented in this paper incorporates several improvements in the 
theoretical calculations and errors compared with previous determinations by OPAL based 
on NLLA+(P(Q;g) QCD. The principal changes are as follows: 

• Improvements have been made in the NLLA theory predictions for the jet broaden- 
ings i?T and -Bw (I’ef. [41] superseding ref. [44]), and for y^^ (ref. [42], superseding 
ref. [45], which was in turn an improvement on ref. [46]). 

• The NLLA calculations for the O-parameter [43] were not available at the time of 
the earliest OPAL publications. 

• The NLLA resummations do not automatically force each event shape distribution 
to vanish at the edge of phase space; missing subleading terms can result in a non¬ 
zero prediction outside the kinematically allowed range of the observable. A remedy 
for this situation involves the substitution 

L = lnfi^ ^ L = \n(--^ + l] . (15) 

\yj \y I/max / 

This method [40] was known at the time of the original OPAL LEP I analysis [2], 
and was investigated as an alternative to the unmodihed NLLA prediction. However, 
this is the hrst time it has been adopted as the standard for as measurements by 
the OPAL Collaboration. 

• The hxed order coefficients A{y) and B{y) are now computed using the EVENT2 
Monte Carlo program [47]. This superseded an earlier program, EVENT [48]. Al¬ 
though both programs are based on the same 0{ag) matrix elements from ref. [38], 
EVENT2 includes an improved algorithm to handle cancellations between real and 
virtual processes, and therefore permits a more precise determination of the coeffi¬ 
cients. 

• A more sophisticated technique for assessing the theoretical errors associated with 
missing higher order terms in the theory has been adopted, based on the extensive 
studies carried out within the LEP QCD working group [39]; see Sect. 5.3 for further 
details. 

The theoretical calculations described above provide predictions of parton-level distri¬ 
butions, i.e. distributions of quarks and gluons. In contrast, the data are corrected to the 
hadron-level, i.e. they correspond to the distributions of the stable particles (including 
photons and both charged and neutral leptons) in the event as explained in Sect. 5.2. 
In order to confront the theory with the hadron-level data, it is necessary to correct the 
theory for the effects of soft fragmentation and hadronization. This was done using large 
samples of events (typically 10^ events) generated using the parton-shower Monte Carlo 
programs PYTHIA (used by default), HERWIG and ARIADNE (used for systematic er¬ 
ror estimates). The analytical theoretical predictions for the cumulative distribution R{y) 
were multiplied by the Monte Carlo prediction for the ratio of R{y) at hadron-level to 
R{y) at parton-leveh. 

^It should be noted that the resummed theoretical calculations apply for massless quarks, while the 
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4.3 QCD predictions for event shape moments 

The moments of the distribution of an event shape observable y are dehned by 


(?/”) 



nl da , 

y -:r^y > 

a dy 


( 16 ) 


where |/max is the kinematically allowed upper limit of the observable. The calculations 
always involve a full integration over phase space, which implies that comparison with data 
always probes all of the available phase space. This is in contrast to QCD predictions for 
distributions; these are commonly only compared with data in restricted regions, where 
the theory is able to describe the data well. Comparisons of QCD predictions for moments 
of event shape distributions with data are thus complementary to tests of the theory using 
the differential distributions. 

The formula for the O^a^) QCD prediction of {y"') is 


iy") = Aa. + 


(17) 


involving the O^a^) coefficients An and 0{a^) coefficients Bn- The values of the coefficients 
An and Bn can be obtained in the same way as described above by running EVENT2. 
The renormalization scale dependence and correction to atot. is implemented in the same 
way as described above for distributions. 

The QCD predictions were transformed from the parton- to the hadron-level by multi¬ 
plying by the ratio of Monte Carlo predictions for the values of moments at the hadron- and 
parton-level. As for the distributions, PYTHIA was used as the standard with HERWIG 
and ARIADNE employed for the estimation of systematic errors. 


5 Analysis methods 

5.1 Selection of events 

The selection of events for this analysis consists of three main stages: the identihcation of 
hadronic event candidates, the removal of events with a large amount of initial-state radi¬ 
ation (ISR) for events at 130 GeV and above, and the removal of four-fermion background 
events for events above 160 GeV, i.e. above the W’''W“ production threshold. 

The selection of hadronic events was based on simple cuts on event multiplicity (to 
remove leptonic hnal states) and on visible energy and longitudinal momentum balance 
(to remove two-photon events). The cuts used at 91 GeV are documented in ref. [49], 
while those used for higher c.m. energies have reoptimized cut parameters, as described in 
ref. [50]. Those parts of the OPAL detector crucial for the present analysis (electromag¬ 
netic calorimeter, jet chamber and trigger system) were required to be fully operational. 

Standard criteria were used to select good tracks and calorimeter energy clusters for 
subsequent analysis. Charged particle tracks were required to have at least 40 hits in the 

quarks in the Monte Carlo models do have masses. No attempt was made to correct for this. In previous 
OPAL papers [2-5], a systematic error was estimated for this effect, and proved to be smaller than the 
other hadronization errors. For consistency with the other LEP experiments, we now neglect this. 
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jet chamber, and at least 50% of the maximum possible number of hits given the polar 
angle of the track. The momentum transverse to the beam axis was required to be at least 
0.15 GeV. Furthermore, the point of closest approach of the track to the collision axis was 
required to be less than 2 cm from the nominal collision point in the x-y plane and less 
than 25 cm in the ; 2 -direction. Energy clusters in the electromagnetic calorimeter were 
required to have energies exceeding 0.10 GeV (0.25 GeV) in the barrel (endcap) region 
of the detector. The number of good charged particle tracks was required to be greater 
than six. After the above cuts the and two-photon background was negligible. 

Furthermore the polar angle of the thrust axis was required to satisfy | cos 6^x1 < 0.9 in 
order that the events be well contained in the detector acceptance. 

At energies signihcantly above Mz, the process of “radiative return” to the Z is a 
common occurrence. In order to study the properties of hadronic events at a well-dehned 
energy scale, it is necessary to eliminate events in which a large amount of energy has been 
lost to ISR. The effective centre-of-mass energy after ISR, was estimated for each 
selected event using the algorithm described in ref. [50]. At c.m. energies of 130 GeV and 
above, we demanded that y/s — y/^ < 10 GeV in order to select a sample of predominantly 
non-radiative events. 

At energies above the W^W“ production threshold, four-fermion events, especially 
those involving qqqq hnal states, become a substantial background. These are reduced 
by using the standard OPAL W'''W“ selection procedure, which is based on a relative 
likelihood method [51]. The same likelihood technique has been applied at each c.m. 
energy studied, with the underlying reference distributions used as inputs to the likelihood 
calculation recomputed for each energy or range of energies. At c.m. energies of 161 GeV 
and above, the qqqq likelihood was required to satisfy /Iqqqq < 0.25 and the qq£z/ likelihood 
was required to satisfy Cqq£u < 0.5. 

After applying the above cuts, the numbers of selected non-radiative qq candidate 
events were as given in Table 1, and were consistent with expectations based on Monte 
Garlo simulated events^. After all cuts, the acceptance for non-radiative qq events (dehned 
for this purpose as those having y/s — y/^ < 1 GeV) ranged from 88.5% at 91 GeV (where 
the loss in acceptance is largely geometrical, arising from the | cos^tI < 0.9 requirement) 
to 76.5% at 207 GeV. The residual four-fermion background was negligible below 161 GeV, 
and otherwise increased from 2.1% at 161 GeV to 6.2% at 207 GeV. 

5.2 Correction procedure 

For each accepted event, the value of each of the event shape observables was computed. 
In order to mitigate the effects of double counting of energy in tracking and calorimetry, a 
standard algorithm was adopted which associated charged particle tracks with calorime¬ 
ter clusters, and subtracted the estimated contribution of the charged particles from the 
calorimeter energy. All selected tracks, and the electromagnetic calorimeter clusters re¬ 
maining after this procedure, were used in the evaluation of event shapes. The event 
shapes were then formed into histograms at each c.m. energy point. In the cases where 

^The numbers of events expected on the basis of Monte Carlo simulations are given in all cases except 
for 91 GeV; to perform an accurate prediction close to the Z peak would require a much more careful 
investigation of the beam energy and luminosity than is required for the present analysis. 
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data at more than one c.m. energy were combined (e.g. 130 and 136 GeV), the data were 
simply summed, and corresponding Monte Carlo samples were created by combining sam¬ 
ples generated at each energy weighted to correspond to the integrated luminosities in 
data. 

The expected number of residual four-fermion background events, was then sub¬ 
tracted from the number of data events, iVj, in each bin, i, of each distribution. The effects 
of detector acceptance and resolution and of residual ISR were then accounted for by a 
simple bin-by-bin correction procedure. For this procedure to be valid, it is necessary 
that the Monte Carlo model give a good description of the data and that the bin size 
be sufficiently large that bin-to-bin migration is reasonably small and symmetric; these 
conditions are sufficiently well satisfied for the present analysis. Two event shape distri¬ 
butions were formed for Monte Carlo simulated qq events; the hrst, at the detector-level, 
treated the Monte Carlo identically to the data, while the second, at the hadron-level, was 
computed using the true momenta of the stable particles in the event^, and was restricted 
to events whose true s' satished ^/s — < 1 GeV. The ratio of the hadron-level to 

the detector-level for each bin, ctj, was used as a correction factor for the data, yielding 
the corrected bin content iV* = ai{Ni — bi). This corrected hadron-level distribution was 
then normalized to unity: Pi = Ni/N , where the sum N = includes any underflow 

and overflow bins. Finally, the differential distribution R'{y) was computed by dividing 
Pi by the bin width. The covariance matrix V for Pi was computed by transforming the 
diagonal Poisson covariance matrix for the uncorrected data iVj: 

Vij = ^2 ^ ^ [Ndik — Nij yNSjk — Njj . (18) 

k ^ ^ k 

By setting o;* = 1 and bi = 0, i.e. iVj = iVj, one may obtain the familiar expression for the 
covariance matrix of a multinomial distribution^. 

The moments were calculated by accumulating sums over all selected events: 

N 

N{y") = J2a • ( 19 ) 

i=l 

where N is the number of events. These sums were corrected by subtracting the back¬ 
ground contribution estimated from simulated events. Then the correction for experi¬ 
mental and ISR effects was performed by multiplying by a detector correction coefficient 
obtained by taking the ratio of the moments at the detector and the hadron-level using 
simulated signal events. The statistical errors of the moments were calculated at the 
detector-level including the effects of Monte Carlo statistics and were subjected to the 
same detector correction. 


"‘For this purpose, all particles having proper lifetimes greater than 3x10 s were regarded as stable. 
^We note therefore that the covariance matrix of a multinomial distribution given by Ni, which has 
sometimes been used in this context, is not identical to the result shown above, and hence not strictly 
correct. 
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5.3 Systematic uncertainties 

Contributions to the systematic uncertainties affecting the corrected hadron-level distri¬ 
butions and moments in data were estimated by repeating the analysis with varied cuts 
or procedures. In each case, the difference in each bin (or for each moment) with respect 
to the standard analysis was taken as a contribution to the systematic error. 

• The containment cut was tightened to | cos6*t| < 0.7. 

• The algorithm to compute s' was replaced by a simpler version in which at most 
one initial state photon was accounted for. 

• The event shapes were computed using all tracks and electromagnetic calorimeter 
clusters. The effects of double counting are then fully taken into account through 
the detector correction procedure. 

• The bin-by-bin corrections, were computed using HERWIG instead of PYTHIA 
as the Monte Carlo hadronization model. 

• The cut on the qqqq four fermion likelihood was tightened to /Iqqqq < 0.1 and 
loosened to /Iqqqq < 0.4; the larger change resulting was taken to be a systematic 
error. 

• The cut on the qqt'z/ four fermion likelihood was tightened to < 0.25 and 

loosened to Cqqei, < 0.75; the larger change resulting was taken to be a systematic 
error. 

• The total four fermion background was varied by ±5%. 

The various contributions above, together with the statistical error on the correction 
factors arising from limited Monte Carlo statistics, were summed in quadrature to form 
the systematic error. None of these sources of systematic error is dominant, but typically 
the larger contributions arise from the use of HERWIG and the use of all tracks and 
clusters. In addition, at high Ec.m. and high y, the variation of the /Iqqqq cut is sometimes 
signihcant. 

In assigning systematic errors to determinations of as, all of the above contributions 
were taken into account and are collectively referred to as the experimental errors. In 
addition, two further sources of systematic error were considered. 

• As explained above, when comparing QCD with the data it is necessary to correct 
for the effects of hadronization. The uncertainty associated with this hadroniza¬ 
tion correction was assessed by using HERWIG 6.2 and ARIADNE 4.11 instead of 
PYTHIA. The larger change in as resulting from these two alternatives was taken 
to define the error. It should be noted that these models have already been tuned 
to similar data to those used here, and hence we adopt this, arguably conservative, 
prescription for assessing the error. 


14 



• The theoretical error, associated with missing higher order terms in the theory, has 
traditionally been assessed by varying the renormalization scale factor, x^, described 
in Sect. 4.2. The predictions of a complete QCD calculation would be independent 
of but a hnite-order calculation such as that used here retains some dependence 
on Xf^. In previous OPAL analyses a range \ < x^ < 2 has been used, and this 
is the procedure we adopt for the analysis of moments. Recently, extensive studies 
have been carried out within the LEP QCD working group [39], which led to a 
more elaborate procedure being proposed, which addresses the question of missing 
higher orders in further ways. We adopt this procedure here for the analysis of 
distributions. A generalization of equation (15) may be considered 


L = 




( 20 ) 


in which p determines how sharply the kinematic cutoff is applied, and x^ acts as 
a scale factor on the event shape. In addition to varying as above, these new 
parameters are varied in the ranges 1 < p < 2 and | < < | (or | < < | 

for An additional matching scheme (i?-matching) is also considered. The 

maximal uncertainty encompassing any of these variations in the theory is taken 
as the systematic error. We take the average of the upper and lower uncertainty 
bands, as dehned in Ref. [39], to dehne the theoretical error. Ref. [39] should be 
consulted for further details. 


6 Results 

6.1 Event shape distributions 

The measured normalized differential cross-sections, R\y), for each of the 14 event shapes 
in each of four energy ranges are given in Tables 2-6.^ The measurements are also shown 
in graphical form in Figures 1-14. In order to clarify the presentation, the data from 161 to 
183 GeV have been combined®, weighted by the numbers of events, and likewise the data 
at 189 GeV and above have been combined in order to provide a high energy set of data 
with reasonably high statistics. These two samples of data, which correspond to mean c.m. 
energies of 177.4 and 197.0 GeV respectively, cover sufficiently small ranges of Ec,m. that 
hard and soft QGD effects should not vary greatly. The horizontal placement of the data 
points follows the prescription in ref. [52], except in a few cases where the 133 GeV and/or 
177 GeV data points have been slightly displaced sideways to avoid overlap. Superimposed 
on Figures 1-14 we show the distributions predicted by the PYTHIA 6.1, HERWIG 6.2 
and ARIADNE 4.11 parton shower models, which in all cases were tuned to other OPAL 

®We note that varies quadratically with the transverse momentum /ct of a radiated gluon, while 
the other observables vary linearly, so by applying this different scaling in the case of we ensure that 
In(fcT) is rescaled by the same amount for all observables. 

^Further details of the data will be made available in the HEPDATA database, 
http://durpdg.dur.ac.uk/HEPDATA/ 

®This combination also has the advantage of reducing any statistical contributions which may be 
present in the systematic error estimates. 
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data recorded during the LEP I running at 91 GeV c.m. energy. The new data appear 
to be well described by all the models. 

In order to make a clearer comparison between data and models, the lower panels 
of Figures 1-14 show the differences between data and each model, divided by the total 
(i.e. statistical and experimental) error on the data in that bin. These ratios are shown 
for 91 GeV and for the combined high energy data sample at 189 GeV and above. The 
sum of squares of these differences would, in the absence of correlations, represent a 
between data and the model. However, since correlations are present, such values 
could be regarded as giving only a rough indication of the agreement between data and 
the models. All three models are seen to describe the high energy data well. Some 
discrepancies are, however, seen in the more precise 91 GeV data. In those observables 
dominated by three-jet prodnction, the largest differences are seen when the observable is 
close to zero, i.e. in the extreme two-jet region. HERWIG is generally apt to give larger 
valnes than the other two models, including some contribution from the extreme three-jet 
regions also. The observables dominated by fonr-jet production (Tmin., Ml, Rn, D and A) 
are less well modelled, and there is a tendency for ARIADNE to give the best description 
of the data. 


6.2 Moments of event shapes 

The measurements of the first hve moments of the event shape observables 1 — T, Mh, 
G, Rt, Rw, ^maj., ^min., S, O, Mi^ and Rn for the four energy ranges are shown in 
Tables 7 and 8.® The same data are shown in Figures 15-18 compared with the same 
Monte Garlo event generators as for the distributions. The lower parts of the fignres show 
again the differences between data and model predictions divided by the total errors. One 
observes that for HERWIG the higher moments generally exhibit larger disagreements. 
This observation is consistent with the distributions where HERWIG showed the most 
signihcant disagreement in the three-jet regions. The PYTHIA and ARIADNE models 
tend to give a better description of the data, the ARIADNE model being somewhat 
closer to the data than PYTHIA. For most observables HERWIG lies above the data at 
91 GeV bnt below at higher energies. The experimental precision of the 91 GeV sample 
is much better than that of the other data samples and thus comparison between data 
and simulation at 91 GeV is more sensitive. 

In order to give an illustration of the sensitivity of the data to QGD effects like the 
running of the strong conpling and the changes in hadronization we compare the hrst 
moments of 1 — T and C measured at 91 and at 197 GeV, see Table 7. The two values of 
(1 — T) are seen to differ by 5.8 standard deviations, treating the experimental errors as 
nncorrelated between the measnrements. Using (C) we observe an 8.6 standard deviation 
difference between the measnrements at 91 and 197 GeV. This shows that onr data are 
indeed sensitive to pertnrbative QGD effects. 


^Further details of the moments will be made available in the HEPDATA database, 
http://durpdg.dur.ac.uk/HEPDATA/ 
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6.3 Determination of CKg 

6.3.1 Determination of Og using event shape distributions 

Our measurement of the strong coupling strength Og is based on hts of QCD predic¬ 
tions to the corrected distributions for 1 — T, Mh, O, -Bt, -Bw and The theoretical 
descriptions of these six observables are among the most complete, allowing the use of 
combined (9(Q;g)-|-NLLA QCD calculations [36,40-43,45]. The main improvements to 
the calculations compared to those used in our previous publications were outlined in 
Sect. 4.2. 

The value of as was estimated by comparing theory with data using a minimum-y^ 
method. In the computation of only the full statistical covariance matrix for the data, 
calculated as explained in Sect. 5.2, was used. Separate hts were performed to each of the 
six observables at each c.m. energy value or range. The ht ranges were the same as those 
used in the previous OPAL publication [5], and were constrained by the requirement that 
the detector and hadronization corrections be reasonably small in the ht region, and that 
both and the htted value of ttg be reasonably stable under small variations in the ht 
range. 

The ht results are summarized in Table 9 for the four c.m. energy points presented 
previously. In addition, further ht results are presented in Table 10 for various other 
groupings of the data in c.m. energy^^. The statistical error on the htted value of as was 
estimated from the variation of by ±1 about its minimum. Systematic uncertainties 
were assessed using the techniques described in Sect. 5.3. For each variant of the analysis, 
the corresponding distribution was htted to determine as, and the diherence with respect 
to the value of as from the default analysis was taken as a systematic error contribution. 
In Figs. 19 and 20 we show the ratio of the data to the htted theory for each of the six 
event shapes at 91 GeV and 197 GeV respectively. Because of normalization, the theory 
predictions are seen to “pivot” about some value of y, which indicates that the data at 
that point have no sensitivity to ag. 

The measurements of as for each observable and c.m. energy are shown in Figure 21. 
We note that the values of cts at 91 GeV are signihcantly higher than at LEP II energies, 
providing evidence for the running of ag. Systematic diherences between the values of 
(Us obtained from diherent observables are seen; for example, (1 — T) and i?T tend to 
give higher than average values of (Ug, whilst i?w tends to give the lowest value. These 
differences may be ascribed to the differing effects of uncomputed higher order terms on 
the various observables; they are often greater than the statistical errors, but are covered 
by the systematic uncertainties. There are also signihcant statistical correlations between 
the values of ctg obtained from the different observables at a given energy, so that the 

^°Note that ref. [16] also determines Ug using there called 02- The small differences between the 
results have been investigated in detail, and are not significant. They may be attributed to differences in 
fit regions, the use of statistically different Monte Carlo samples, and the adoption of slightly different 
strategies for the assessment of theoretical errors. 

^^These energy values and ranges are those used by the LEP QCD working group, namely 161 GeV, 
172 GeV, 183 GeV, 189-192 GeV, 196-202 GeV and > 202 GeV. In the case of the 196-202 GeV point, 
the value of Ug has been run from the mean c.m. energy of 198.6 GeV to a nominal value of 200 GeV 
using the expected QGD behaviour — a correction of 0.0001. 
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scatter of the points at energies where the statistics are low (e.g. 161 GeV) tends to be 
smaller than one would expect if the statistical errors were uncorrelated. 

It is useful to combine the measurements of as from different observables and/or c.m. 
energy points, in order to exploit the data fully. This problem has been the subject of 
extensive study within the LEP QCD working group [53], and we adopt their procedure 
here. In brief, the method is as follows. The as measurements to be combined are hrst 
evolved to a common scale, Q = Mz, assuming the validity of QCD. Their combination is 
then performed using a weighted mean method, such as to minimize the between the 
combined value and the measurements. So, if the measured values evolved to Q = Mz 
are denoted with covariance matrix V, the combined value, as{Mz), is given by 

as{Mz) = ^Wiai where Wi = ■ (21) 

i j j,k 

The combined values may then be evolved back to the original scale if required. The 
difficulty resides in making a reliable estimate of V in the presence of dominant and highly 
correlated systematic errors; small uncertainties in the estimation of these correlations can 
easily cause undesirable features such as negative weights. For this reason, only statistical 
errors (estimated using data-sized subsamples of Monte Carlo events to assess correlations 
between different observables in a given data sample) and experimental systematic errors 
(for which the correlations are estimated using the “minimum overlap” assumption^^) 
were allowed to contribute to the off-diagonal elements in V when computing the weights, 
while all error contributions were included in the diagonal terms. The hadronization and 
theoretical uncertainties were computed by combining the as values obtained with the 
alternative hadronization models, and from the upper and lower theoretical errors, using 
the same weights. 

If the full covariance matrix were used to compute the weights Wi, then this choice 
of weights would have the effect of minimising the total error. However, the correlations 
between the systematic errors cannot be completely reliably estimated (as manifested by 
the occurrence of negative weights). The modified procedure adopted here will not in 
general minimize the error, and indeed the error on the weighted average may be greater 
than one of the inputs; this actually arises in our data for the case of Despite this, 
we consider the combined value to provide the safest estimate of ag, since it cannot be 
guaranteed that the relatively small theoretical error for is not fortuitous. 

In Table 11 we give the values of as for each observable, evolved to a common scale 
Mz, combined over all c.m. energies. These results are also summarized in Figure 22. 
This shows clearly that the measurements from the different observables are far from 
compatible when only statistical errors are considered, but are consistent with a common 
mean when the systematic errors are included. The results of combining the as values for 
the six observables are given in the rightmost columns of Tables 9 and 10. In addition, 
the relative weight, Wi, carried by each observable is given, from which we see that 
generally carries the greatest weight (because it has the smallest theoretical uncertainty). 
These results are plotted in Figure 23 where we show the values at each energy evolved to 

^^The minimum overlap ansatz involves taking the covariance between any pair of systematic error 
contributions to be equal to the smaller of the two variances. 
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a common scale, as{Mz), and in Figure 24 where we show Os as a function of the energy 
scale Q = -Fc.m.- These plots show that the variation of ttg with -Fc.m. is consistent with 
the running predicted by QCD, and is incompatible with a constant value of ag. For 
example, the two most precise values of ag, at 91 and 197 GeV, differ by more than three 
standard deviations (applying the minimum overlap ansatz to the systematic errors). 

The measurement of cts based on the 91 GeV data is 

ag{Mz) = 0.1192 ± 0.0002 (stat.) ± 0.0008 (expt.) ± 0.0015 (hadr.) ± 0.0047 (theo.) 

while the result for ag{Mz) combining all higher energy data at Q > Mz is 

ag{Mz) = 0.1189 ± 0.0011 (stat.) ± 0.0015 (expt.) ± 0.0008 (hadr.) ± 0.0040 (theo.) . 

The consistency between these measurements, which are equal within the statistical errors, 
again shows that the data are compatible with the running predicted by QGD, since this 
running was assumed in evolving the high energy data to Q = Mz- We note that the high 
energy data have signihcantly smaller theoretical and hadronization errors, and therefore 
complement the statistically superior 91 GeV data. The value of as{Mz) obtained from 
all observables and all energies combined is^^: 

ag{Mz) = 0.1191 ± 0.0005 (stat.) ± 0.0010 (expt.) ± 0.0011 (hadr.) ± 0.0044 (theo.) 

= 0.1191 ± 0.0005 (stat.) ± 0.0046 (syst.) 

= 0.1191 ±0.0047 (total) 


6.3.2 Determination of otg using event shape moments 

The strong coupling ag{Mz) has also been determined from the measured moments using 
the 0(ag) QGD predictions explained in Sect. 4.3. This is the first such study to be 
published by the OPAL Gollaboration. The hrst hve moments of the observables 1 — T, 
Mh, Bt, By/, C and were studied, i.e. the same observables as used for the deter¬ 
mination of Os from distributions. One might anticipate a priori that this would yield 
a less precise determination of ag than the differential distributions, because the theory 
lacks resummation of large logarithms, and the moments include regions of phase space 
where hadronization effects are large. Nevertheless, the comparison of ag determined in 
this way with that obtained from the distributions should provide an illuminating test of 
the adequacy of QGD in this area. 

The fits proceeded by comparing the data at the four combined energy points for 
a given moment of an observable with the theory prediction. The running of ag was 
implemented in the £t in two-loop precision using the formula given in ref. [54]. A value 
of was calculated using the statistical errors of the data and minimized to extract a 
value of ag{Mz)- The hts were repeated for each systematic variation of the analysis. The 
statistical error was found as above by variation of by ±1 about the minimum. 

The £t results are shown in Table 12 and in Figure 25. The fit to (Mh) was not stable 
and therefore no results are shown. We obtain values of y^/d.o.f. of 0(1); the fitted QGD 

^^The x^/d.o.f. for this overall fit is 10.5/23. The value is significantly smaller than unity because of 
the necessity to neglect correlations in the combination. 
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predictions including the running of ag are thus consistent with our data. However, we 
find that for ((1 —T)”), (C”) and (B^) the htted values of as{Mz) increase with the 
order n of the moment used. This effect is not observed for (B^), and (Mg), 

n = 2,..., 5. In order to investigate the origin of this behaviour we show in Figure 26 
the ratio K = Bn/An of NLO and LO coefficients for the six observables used in our hts. 
There is a clear correlation between the increasing values of as{Mz) with moment n and 
increasing values of K with n for ((1 — T)'^), (C") and {B^). The other observables (B^), 
((y^)”) and (Mg), n = 2,..., 5, have fairly constant values of K and correspondingly 
stable results for as{Mz). We also note that (Mh) has a large and negative value of K 
which is the cause of the unstable hts^^. 

In order to combine the individual determinations of ag{Mz) we used the same weighted 
average method as in Sect. 6.3.1. We considered only those results for which the £t was 
successful and for which the NLO term in equation (17) is less than half the LO term (i.e. 
\Kag/2'K\ < 0.5 or \K\ < 25), namely (1 - T), (O), (Ht), (5(^) and ((i/^g)’"), n = 1,..., 5 
and (Mg), n = 2,..., 5; i.e. results from 17 observables in total. The statistical cor¬ 
relations between the 17 results were determined using Monte Carlo simulation at the 
hadron-level. The experimental errors were included in the covariance matrix using the 
minimum-overlap assumption while the hadronization and renormalization scale uncer¬ 
tainties were only added to the diagonal of the covariance matrix. The combination was 
repeated using the same weights with the HERWIG and ARIADNE hadronization cor¬ 
rections and for = 0.5 and Xfj, = 2.0 for the calculation of the hadronization theory 
systematic errors. The weights were generally (9(10%) for each observable; the largest 
weight was 23% for (Ht) while the only weights below 3% were those for ((y^)^) (Mg) 
and (Mg). The result is 

ag(Mz) = 0.1223 ± 0.0005 (stat.) ± 0.0014 (expt.) ± 0.0016 (hadr.)gg“®R(theo.) 

= 0.1223 ± 0.0005 (stat.) ± 0.0058 (syst.) 

= 0.1223 ± 0.0059 (total) 

in good agreement with the result from distributions presented in Sect. 6.3.1. 

The experimental, hadronization and theory uncertainties are somewhat larger than 
for the distributions. In an analysis using moments the complete available phase space is 
sampled including regions which are more difficult to measure experimentally and which 
are less reliably modelled by the hadronization models. Also, the NLO QCD prediction is 
less complete than the matched (9(Q;g)-|-NLLA prediction available for the distributions 
which we studied. It is nevertheless a remarkable success of QCD together with the 
corresponding hadronization models that the NLO theory is able to describe successfully 
some observables based on the complete available phase space. 

7 Summary and conclusions 

In this paper we have presented measurements of the event shapes for hadronic events 
produced at LEP at centre-of-mass energies between 91 and 209 GeV. Both differential 

^■^By reference to equation (17) we can see that there is no real solution for as if Bn < -Al/Ay-). 
This is the case for (Mh), which accounts for the failure of the hts. There is one other case where the 
coefficient Bn is negative, namely (Uw), but it is not sufficiently negative to preclude a solution for a^. 
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distributions and moments have been determined. 

The predictions of the PYTHIA, HERWIG and ARIADNE Monte Carlo models are 
found to be in general agreement with the measured distributions and moments. However, 
some discrepancies are noted in the 91 GeV data where the statistical errors are smallest. 
The main differences between models and data occur in the extreme two-jet region, and 
for observables sensitive to the production of four or more jets. In general, ARIADNE 
provides the best description of the data and HERWIG the least good. 

From a £t of (P(Q;g)-|-NLLA QCD predictions to the distributions of six event shape 
observables, we have determined the strong coupling parameter ag. The variation of ag 
with energy scale over the range 91 to 209 GeV is found to be in accordance with the 
expectations of QCD. For example, the measurement of ag based on the 91 GeV data is 

ag{Mz) = 0.1192 ± 0.0002 (stat.) ± 0.0008 (expt.) ± 0.0015 (hadr.) ± 0.0047 (theo.) . 

Assuming the validity of QCD, the higher energy measurements can all be evolved to a 
common scale Q = Mz and combined, yielding the following result for ag{Mz) combining 
all higher energy data at Q > Mz 

ag{Mz) = 0.1189 ± 0.0011 (stat.) ± 0.0015 (expt.) ± 0.0008 (hadr.) ± 0.0040 (theo.) . 

Combining all data at all energy scales, the value of ag{Mz) is determined to be 

ag{Mz) = 0.1191 ± 0.0005 (stat.) ± 0.0010 (expt.) ± 0.0011 (hadr.) ± 0.0044 (theo.) , 

in good agreement with the world average quoted in ref. [54]. The results for ag{Mz) 
derived from different event shapes are consistent within errors. 

Values of ag{Mz) have also been derived from the energy dependence of event shape 
moments, using O(a^) QCD. Although less complete than the (P(Q;g)-|-NLLA QCD pre¬ 
dictions used for the distributions, these calculations prove to give a consistent description 
of many, though not all, of the moments. The combined value obtained from the moments 
was ag{Mz) = 0.1223 ± 0.0005 (stat.) ± 0.0014 (expt.) ± 0.0016 (hadr.))'(QQQ 3 g(theo.) , con¬ 
sistent with that derived from the distributions. However, because the value obtained from 
the distributions is based on the more complete theory (and has a smaller overall error) 
we consider that to be the most reliable estimate from the present data. 
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Year 

Range of 
(GeV) 

Mean ^/s 
(GeV) 

Integrated 

Inminosity 

(pb-i) 

Nnmber of 
selected 

events 

Expected 

nnmber 

1996-2000 

91.0-91.5 

91.3 

14.7 

395695 

— 

1995, 1997 

129.9- 136.3 

133.1 

11.26 

630 

698 


161.2-161.6 

161.3 

10.06 

281 

275 

HH 

170.2- 172.5 

172.1 

10.38 

218 

232 

■H 

180.8- 184.2 

182.7 

57.72 

1077 

1084 

1998 

188.3- 189.1 

188.6 

185.2 

3086 

3130 

1999 

191.4-192.1 

191.6 

29.53 

514 

473 

1999 

195.4- 196.1 

195.5 

76.67 

1137 

1161 

1999, 2000 

199.1 - 200.2 

199.5 

79.27 

1090 

1131 

1999, 2000 

201.3- 202.1 

201.6 

37.75 

519 

527 

2000 

202.5- 205.5 

204.9 

82.01 

1130 

1090 

2000 

205.5- 208.9 

206.6 

138.8 

1717 

1804 


Table 1: The OPAL data samples used for the present analysis. The horizontal lines 
separate the data into the four energy ranges used for presentation purposes. 
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1 - 

T 

R'il 

- 

T ) at 91 GeV 

R '{1 

- r) at 133 GeV 

R'(l 

- T ) at 177 GeV 

R '( l - 

- T ) at 197 GeV 

0.00 - 

0.01 

1.276 

± 

0.019 

± 

0.040 

4.37 

± 

0.82 

± 

0.55 

8.76 

± 0.74 

± 0.64 

9.95 

± 0.34 

± 0.27 

0.01 - 

0.02 

12.25 

± 

0.05 

± 

0.41 

20.4 

± 

1.6 

± 

3.8 

22.65 

± 1.11 

± 1.23 

22.44 

± 0.47 

± 0.40 

0.02 - 

0.03 

18.38 

± 

0.06 

± 

0.27 

20.4 

± 

1.5 

± 

1.7 

16.25 

± 0.92 

± 0.76 

14.92 

± 0.37 

± 0.56 

0.03 - 

0.04 

13.87 

± 

0.06 

± 

0.14 

10.5 

± 

1.2 

± 

0.9 

9.84 

± 0.73 

± 0.98 

9.57 

± 0.31 

± 0.41 

0 . 04 - 

0.05 

9.83 

± 

0.05 

± 

0.16 

6.7 

± 

1.0 

± 

0.4 

6.85 

± 0.65 

± 0.87 

7.23 

± 0.27 

± 0.20 

0.05 - 

0.07 

6.502 

± 

0.026 

± 

0.088 

4.70 

± 

0.58 

± 

0.35 

4.96 

± 0.38 

± 0.20 

5.24 

± 0.16 

± 0.15 

0 . 07 - 

0.09 

4.127 

± 

0.022 

± 

0.041 

3.64 

± 

0.51 

± 

0.29 

3.18 

± 0.33 

± 0.27 

2.87 

± 0.14 

± 0.15 

0.09 - 

0.12 

2.646 

± 

0.014 

± 

0.058 

1.68 

± 

0.32 

± 

0.44 

2.05 

± 0.23 

± 0.19 

2.28 

± 0.09 

± 0.07 

0.12 - 

0.15 

1.709 

± 

0.012 

± 

0.079 

1.36 

± 

0.24 

± 

0.32 

1.28 

± 0.19 

± 0.21 

1.40 

± 0.08 

± 0.10 

0.15 - 

0.22 

0.910 

± 

0.005 

± 

0.019 

1.09 

± 

0.12 

± 

0.18 

0.82 

± 0.10 

± 0.10 

0.787 ± 0.047 ± 0.057 

0.22 - 

0.30 

0.3705 ± 

0.0032 ± 

0.0092 

0.467 ± 0.065 ± 

0.066 

0.342 ± 0.075 ± 0.088 

0.382 ± 0.046 ± 0.065 

Mh 

R ' (Mh) at 91 GeV 

R ' (Mh) at 133 GeV 

R ' (Mh) at 177 GeV 

R ' (Mh) at 197 GeV 

0.060 - 

0.075 

0.119 

± 

0.003 

± 

0.010 

0.38 

± 

0.16 

± 

0.13 

1.35 

± 0.18 

± 0.44 

1.93 

± 0.09 

± 0.17 

0.075 - 

0.090 

0.55 

± 

0.01 

± 

0.13 

1.56 

± 

0.35 

± 

0.36 

4.33 

± 0.36 

± 0.58 

4.62 

± 0.17 

± 0.19 

0.090 - 

0.110 

2.16 

± 

0.01 

± 

0.39 

4.69 

± 

0.53 

± 

1.25 

6.21 

± 0.45 

± 0.59 

6.55 

± 0.19 

± 0.23 

0.110 - 

0.140 

5.72 

± 

0.02 

± 

0.05 

8.18 

± 

0.54 

± 

1.05 

6.06 

± 0.34 

± 0.43 

5.82 

± 0.14 

± 0.30 

0.140 - 

0.170 

6.66 

± 

0.02 

± 

0.36 

4.99 

± 

0.48 

± 

0.82 

5.06 

± 0.30 

± 0.26 

4.66 

± 0.12 

± 0.16 

0.170 - 

0.200 

4.88 

± 

0.02 

± 

0.15 

3.92 

± 

0.43 

± 

0.28 

3.47 

± 0.26 

± 0.43 

3.53 

± 0.11 

± 0.16 

0.200 - 

0.250 

3.29 

± 

0.01 

± 

0.12 

2.47 

± 

0.27 

± 

0.30 

2.54 

± 0.18 

± 0.36 

2.37 

± 0.07 

± 0.12 

0.250 - 

0.300 

2.107 

± 

0.010 

± 

0.054 

1.41 

± 

0.23 

± 

0.22 

1.54 

± 0.16 

± 0.28 

1.60 

± 0.07 

± 0.09 

0.300 - 

0.350 

1.352 

± 

0.008 

± 

0.042 

1.35 

± 

0.18 

± 

0.23 

1.25 

± 0.14 

± 0.18 

1.22 

± 0.06 

± 0.12 

0.350 - 

0.450 

0.703 

± 

0.004 

± 

0.031 

0.758 ± 

0.088 ± 

0.075 

0.647 ± 0.079 ± 0.085 

0.641 

± 0.035 

± 0.031 

0.450 - 

0.600 

0.1372 ± 

0.0016 ± 

0.0036 

0.162 ± 

0.033 ± 

0.016 

0.114 ± 0.032 ± 0.042 

0.149 

± 0.019 

± 0.039 

C 

R ! ( C ) at 91 GeV 

R ' ( C ) at 133 GeV 

R ' ( C ) at 177 GeV 

R ' ( C ) at 197 GeV 

0.00 - 

0.05 

0.2186 ± 

0.0038 ± 

0.0060 

0.83 

± 

0.17 

± 

0.11 

1.97 

± 0.16 

± 0.12 

2.079 ± 0.071 ± 0.053 

0.05 - 

0.08 

2.061 

± 

0.013 

± 

0.108 

4.07 

± 

0.45 

± 

1.17 

5.04 

± 0.32 

± 0.31 

5.225 ± 0.136 ± 0.159 

0.08 - 

0.11 

4.037 

± 

0.018 

± 

0.052 

5.70 

± 

0.46 

± 

0.51 

4.32 

± 0.28 

± 0.43 

3.797 ± 0.113 ± 0.235 

0.11 - 

0.14 

4.152 

± 

0.018 

± 

0.047 

3.86 

± 

0.39 

± 

0.41 

3.26 

± 0.24 

± 0.14 

3.076 ± 0.099 ± 0.052 

0.14 - 

0.18 

3.225 

± 

0.013 

± 

0.044 

2.55 

± 

0.29 

± 

0.14 

2.30 

± 0.18 

± 0.22 

2.382 ± 0.075 ± 0.045 

0.18 - 

0.22 

2.421 

± 

0.012 

± 

0.060 

1.65 

± 

0.25 

± 

0.19 

1.73 

± 0.16 

± 0.11 

1.757 ± 0.067 ± 0.026 

0.22 - 

0.30 

1.705 

± 

0.007 

± 

0.020 

1.22 

± 

0.15 

± 

0.07 

1.317 ± 0.097 ± 0.085 

1.330 ± 0.041 ± 0.033 

0.30 - 

0.40 

1.112 

± 

0.005 

± 

0.012 

0.93 

± 

0.11 

± 

0.06 

0.853 ± 0.076 ± 0.067 

0.850 ± 0.032 ± 0.028 

0.40 - 

0.50 

0.747 

± 

0.004 

± 

0.017 

0.43 

± 

0.09 

± 

0.12 

0.573 ± 0.067 ± 0.062 

0.623 ± 0.028 ± 0.017 

0.50 - 

0.60 

0.535 

± 

0.004 

± 

0.024 

0.73 

± 

0.08 

± 

0.15 

0.488 ± 0.063 ± 0.072 

0.453 ± 0.027 ± 0.026 

0.60 - 

0.75 

0.3693 ± 

0.0024 ± 

0.0082 

0.334 ± 

0.051 ± 

0.061 

0.318 ± 0.048 ± 0.055 

0.341 ± 0.025 ± 0.047 

0.75 - 

1.00 

0.0982 ± 

0.0009 ± 

0.0023 

0.101 ± 

0.018 ± 

0.021 

0.069 ± 0.029 ± 0.034 

0.070 ± 0.021 ± 0.045 


Table 2: Distributions for the thrust (1 — T), heavy jet mass (Mh) and C-parameter measured by OPAL at centre-of-mass 
energies ^/s = 91, 130-136, 161-183, and 189-209 GeV. The hrst uncertainty is statistical, while the second is systematic. 





R ' ( B ' y ) at 91 GeV 

R ' ( Bt ) at 133 GeV 

R ' ( Bt ) at 177 GeV 

B '( Bt 

) at 197 GeV 

0.000 - 0.030 

0.117 

± 0.005 

± 0.012 

0.58 

± 

0.26 

± 

0.16 

1.78 

± 

0.23 

± 

0.18 

2.40 

± 

0.12 

± 

0.07 

0.030 - 0.040 

2.50 

± 0.03 

± 0.11 

10.2 

± 

1.1 

± 

1.3 

12.1 

± 

0.9 

± 

1.2 

11.23 

± 

0.41 

± 

0.44 

0.040 - 0.050 

6.95 

± 0.04 

± 0.13 

9.6 

± 

1.3 

± 

1.9 

11.5 

± 

0.8 

± 

1.2 

10.99 

± 

0.38 

± 

0.49 

0.050 - 0.060 

9.806 

± 0.050 

± 0.079 

11.9 

± 

1.2 

± 

0.9 

11.01 

± 

0.75 

± 

0.79 

9.67 

± 

0.34 

± 

0.56 

0.060 - 0.075 

10.744 

± 0.040 

± 0.120 

10.06 

± 

0.86 

± 

0.52 

8.29 

± 

0.55 

± 

0.43 

7.74 

± 

0.25 

± 

0.23 

0.075 - 0.090 

8.777 

± 0.035 

± 0.211 

6.31 

± 

0.76 

± 

0.36 

5.71 

± 

0.48 

± 

0.36 

6.14 

± 

0.22 

± 

0.35 

0.090 - 0.110 

6.597 

± 0.026 

± 0.065 

4.35 

± 

0.59 

± 

0.82 

4.38 

± 

0.37 

± 

0.38 

4.85 

± 

0.18 

± 

0.16 

0.110 - 0.130 

4.829 

± 0.023 

± 0.074 

4.36 

± 

0.52 

± 

0.41 

3.64 

± 

0.34 

± 

0.36 

3.43 

± 

0.15 

± 

0.14 

0 . 130 - 0.160 

3.386 

± 0.016 

± 0.050 

2.43 

± 

0.36 

± 

0.25 

2.78 

± 

0.24 

± 

0.30 

2.60 

± 

0.11 

± 

0.08 

0.160 - 0.200 

2.130 

± 0.011 

± 0.047 

1.71 

± 

0.24 

± 

0.22 

1.50 

± 

0.18 

± 

0.17 

1.747 ± 

0.081 ± 

0.063 

0.200 - 0.250 

1.186 

± 0.007 

± 0.025 

1.25 

± 

0.16 

± 

0.30 

1.13 

± 

0.15 

± 

0.13 

1.008 ± 

0.069 ± 

0.085 

0.250 - 0.300 

0.565 

± 0.005 

± 0.015 

0.72 

± 

0.10 

± 

0.09 

0.42 

± 

0.12 

± 

0.11 

0.443 ± 

0.082 ± 

0.073 

0.300 - 0.350 

0.1557 ± 0.0024 ± 0.0049 

0.128 ± 

0.042 ± 

0.059 

0.05 

± 

0.09 

± 

0.13 

0.096 ± 

0.082 ± 

0.114 

Bw 

R ' at 91 GeV 

R ^ (^ w ) ^t 133 GeV 


■) at 177 GeV 


) at 197 GeV 

0.000 - 0.020 

0.557 

± 0.011 

± 0.021 

3.15 

± 

0.53 

± 

0.83 

5.29 

± 

0.45 

± 

0.23 

5.93 

± 

0.19 

± 

0.19 

0.020 - 0.030 

10.38 

± 0.05 

± 0.15 

16.4 

± 

1.5 

± 

2.5 

16.76 

± 

0.94 

± 

1.33 

15.82 

± 

0.38 

± 

0.90 

0.030 - 0.040 

16.54 

± 0.06 

± 0.11 

14.6 

± 

1.3 

± 

0.6 

13.23 

± 

0.83 

± 

0.74 

12.71 

± 

0.34 

± 

0.26 

0.040 - 0.050 

13.32 

± 0.05 

± 0.66 

10.7 

± 

1.1 

± 

0.7 

11.09 

± 

0.75 

± 

0.79 

10.13 

± 

0.31 

± 

0.37 

0.050 - 0.065 

9.82 

± 0.04 

± 0.13 

9.47 

± 

0.84 

± 

0.84 

7.70 

± 

0.54 

± 

0.51 

7.79 

± 

0.22 

± 

0.18 

0.065 - 0.080 

7.17 

± 0.03 

± 0.15 

5.44 

± 

0.72 

± 

0.40 

5.05 

± 

0.47 

± 

0.21 

5.55 

± 

0.20 

± 

0.17 

0.080 - 0.100 

5.061 

± 0.023 

± 0.065 

3.96 

± 

0.54 

± 

0.45 

4.24 

± 

0.37 

± 

0.33 

4.19 

± 

0.16 

± 

0.19 

0.100 - 0.150 

2.845 

± 0.011 

± 0.066 

2.44 

± 

0.27 

± 

0.24 

2.47 

± 

0.19 

± 

0.12 

2.442 ± 

0.079 ± 

0.096 

0.150 - 0.200 

1.238 

± 0.008 

± 0.042 

1.19 

± 

0.16 

± 

0.25 

1.11 

± 

0.15 

± 

0.15 

1.146 ± 

0.066 ± 

0.062 

0.200 - 0.250 

0.465 

± 0.005 

± 0.014 

0.63 

± 

0.11 

± 

0.07 

0.47 

± 

0.11 

± 

0.23 

0.506 ± 

0.057 ± 

0.087 

0.250 - 0.300 

0.0625 ± 0.0017 ± 0.0029 

0.077 ± 

0.033 ± 

0.057 

0.10 

± 

0.05 

± 

0.10 

0.123 ± 

0.031 ± 

0.095 

2/?3 

R ' {yf^) at 91 GeV 

R'{y?3 

) at 133 GeV 

R'{y?3 

) at 177 GeV 

R'iy?3 

) at 197 GeV 

0.00030 - 0.00075 

146.5 

± 1.0 

± 2.6 

299 

± 

30 

± 

41 

323 

± 

20 

± 

19 

334 

± 

8 

± 

17 

0.00075 - 0.00130 

183.1 

± 0.9 

± 5.5 

201 

± 

22 

± 

21 

223 

± 

14 

± 

36 

168.6 

± 

5.5 

± 

5.2 

0.00130 - 0.00230 

141.4 

± 0.6 

± 3.1 

116 

± 

12 

± 

22 

93 

± 

7 

± 

16 

101.5 

± 

3.0 

± 

4.8 

0.00230 - 0.00400 

81.8 

± 0.3 

± 1.3 

62.6 

± 

6.7 

± 

8.4 

53.1 

± 

4.2 

± 

7.1 

58.9 

± 

1.8 

± 

2.6 

0.00400 - 0.00700 

39.6 

± 0.2 

± 1.2 

33.0 

± 

3.9 

± 

5.1 

32.8 

± 

2.6 

± 

5.1 

30.2 

± 

1.0 

± 

0.7 

0.00700 - 0.01200 

19.90 

± 0.09 

± 0.22 

19.9 

± 

2.2 

± 

3.1 

16.7 

± 

1.5 

± 

1.4 

17.59 

± 

0.62 

± 

0.30 

0.01200 - 0.02250 

9.73 

± 0.05 

± 0.20 

6.9 

± 

1.1 

± 

0.4 

8.67 

± 

0.76 

± 

0.92 

8.96 

± 

0.32 

± 

0.44 

0.02250 - 0.04000 

4.56 

± 0.02 

± 0.20 

3.72 

± 

0.60 

± 

0.50 

4.54 

± 

0.41 

± 

0.34 

3.90 

± 

0.18 

± 

0.28 

0.04000 - 0.07000 

2.105 

± 0.013 

± 0.098 

2.30 

± 

0.29 

± 

0.35 

1.94 

± 

0.23 

± 

0.27 

2.08 

± 

0.10 

± 

0.18 

0.07000 - 0.13000 

0.824 

± 0.006 

± 0.024 

0.93 

± 

0.12 

± 

0.25 

0.75 

± 

0.11 

± 

0.09 

0.773 ± 

0.049 ± 

0.033 

0.13000 - 0.23500 

0.266 

± 0.003 

± 0.014 

0.315 ± 

0.058 ± 

0.069 

0.243 ± 

0.054 ± 

0.059 

0.252 ± 

0.029 ± 

0.055 

0.23500 - 0.33333 

0.0506 ± 0.0014 ± 0.0033 

0.034 ± 

0.026 ± 

0.018 

0.023 ± 

0.034 ± 

0.039 

0.056 ± 

0.021 ± 

0.041 


Table 3: Distributions for the total jet broadening (-Bt) wide jet broadening (-Bw) and Durham two- to three-jet transition pa¬ 
rameter (i/^) measured by OPAL at centre-of-mass energies y/s = 91, 130-136, 161-183, and 189-209 GeV. The hrst uncertainty 
is statistical, while the second is systematic. 




T ■ 

R' (Tnaj.) at 91 GeV 

R' (Tnaj.) at 133 GeV 

R' (Tnaj. 

) at 177 GeV 

R' (Tnaj.) at 197 GeV 

0.00 - 0.04 

0.0363 ± 0.0028 ± 0.0073 

0.27 

± 0.15 

± 0.18 

0.94 

± 

0.15 

± 

0.13 

1.185 ± 0.069 ± 0.025 

0.04 - 0.08 

3.373 

± 0.015 

± 0.041 

6.00 

± 0.44 

± 0.60 

6.99 

± 

0.29 

± 

0.29 

6.70 

± 0.12 

± 0.10 

0.08 - 0.12 

6.499 

± 0.018 

± 0.124 

5.93 

± 0.39 

± 0.16 

5.18 

± 

0.25 

± 

0.15 

4.81 

± 0.10 

± 0.15 

0.12 - 0.16 

4.415 

± 0.015 

± 0.065 

3.32 

± 0.33 

± 0.40 

3.27 

± 

0.21 

± 

0.17 

3.34 

± 0.09 

± 0.10 

0.16 - 0.22 

2.753 

± 0.010 

± 0.046 

2.29 

± 0.23 

± 0.18 

2.03 

± 

0.15 

± 

0.17 

2.124 ± 0.062 ± 0.066 

0.22 - 0.30 

1.579 

± 0.007 

± 0.027 

1.12 

± 0.15 

± 0.16 

1.25 

± 

0.11 

± 

0.08 

1.327 ± 0.044 ± 0.051 

0.30 - 0.40 

0.825 

± 0.004 

± 0.023 

0.94 

± 0.10 

± 0.11 

0.74 

± 

0.08 

± 

0.10 

0.722 ± 0.034 ± 0.047 

0.40 - 0.50 

0.3864 ± 0.0030 ± 0.0082 

0.363 ± 0.064 ± 0.059 

0.361 ± 

0.065 ± 

0.059 

0.381 ± 0.033 ± 0.065 

0.50 - 0.60 

0.1348 ± 0.0017 ± 0.0029 

0.220 ± 0.038 ± 0.046 

0.125 ± 

0.046 ± 

0.055 

0.129 ± 0.031 ± 0.046 

Tnin. 

R' (Tnin.) at 91 GeV 

R' (Tnin.) at 133 GeV 

R' (Tnin. 

) at 177 GeV 

R' (Tnin.) at 197 GeV 

0.00 - 0.02 

0.0220 ± 0.0028 ± 0.0029 

0.07 

± 0.08 

± 0.19 

0.154 ± 

0.096 ± 

0.050 

0.336 ± 0.049 ± 0.062 

0.02 - 0.04 

1.388 

± 0.015 

± 0.040 

5.13 

± 0.62 

± 0.44 

9.90 

± 

0.55 

± 

0.65 

9.75 

± 0.28 

± 0.39 

0.04 - 0.06 

8.151 

± 0.031 

± 0.078 

13.9 

± 0.9 

± 1.8 

14.77 

± 

0.66 

± 

0.82 

15.07 

± 0.34 

± 0.50 

0.06 - 0.08 

12.415 

± 0.036 

± 0.101 

12.1 

± 0.8 

± 1.7 

9.81 

± 

0.54 

± 

0.65 

9.39 

± 0.24 

± 0.25 

0.08 - 0.10 

10.342 

± 0.032 

± 0.065 

6.9 

± 0.7 

± 1.1 

5.38 

± 

0.40 

± 

0.34 

5.72 

± 0.17 

± 0.20 

0.10 - 0.12 

6.852 

± 0.027 

± 0.066 

4.69 

± 0.49 

± 0.62 

3.50 

± 

0.32 

± 

0.33 

2.90 

± 0.13 

± 0.22 

0.12 - 0.14 

4.186 

± 0.021 

± 0.043 

2.75 

± 0.38 

± 0.69 

2.36 

± 

0.27 

± 

0.35 

2.00 

± 0.12 

± 0.18 

0.14 - 0.16 

2.423 

± 0.016 

± 0.054 

1.59 

± 0.31 

± 0.38 

1.45 

± 

0.23 

± 

0.26 

1.26 

± 0.11 

± 0.12 

0.16 - 0.20 

1.255 

± 0.008 

± 0.024 

0.94 

± 0.16 

± 0.35 

0.55 

± 

0.14 

± 

0.17 

0.69 

± 0.07 

± 0.06 

0.20 - 0.24 

0.499 

± 0.005 

± 0.017 

0.217 ± 0.095 ± 0.049 

0.05 

± 

0.12 

± 

0.12 

0.45 

± 0.07 

± 0.22 

0.24 - 0.30 

0.1733 ± 0.0023 ± 0.0047 

0.145 ± 0.046 ± 0.029 

0.22 

± 

0.11 

± 

0.18 

0.01 

± 0.07 

± 0.10 

A 

R ' { A ) at 91 GeV 

R ' { A ) at 133 GeV 

R ' (^) at 177 GeV 

R ' ( A ) at 197 GeV 

0.000 - 0.005 

76.9 

± 0.2 

± 1.2 

115.0 

± 3.8 

± 7.4 

126 

± 

4 

± 

10 

135.3 

± 2.5 

± 5.4 

0.005 - 0.010 

54.68 

± 0.14 

± 0.53 

38.3 

± 3.1 

± 2.7 

31.2 

± 

2.2 

± 

2.7 

30.0 

± 0.9 

± 2.3 

0.010 - 0.015 

25.52 

± 0.11 

± 0.43 

16.8 

± 2.0 

± 1.9 

13.0 

± 

1.4 

± 

0.6 

11.8 

± 0.6 

± 1.2 

0.015 - 0.025 

10.90 

± 0.05 

± 0.15 

7.7 

± 0.9 

± 1.7 

6.11 

± 

0.68 

± 

0.68 

4.81 

± 0.31 

± 0.38 

0.025 - 0.040 

3.686 

± 0.022 

± 0.065 

2.60 

± 0.47 

± 0.24 

1.60 

± 

0.40 

± 

0.34 

1.97 

± 0.19 

± 0.33 

0.040 - 0.070 

1.111 

± 0.008 

± 0.021 

0.76 

± 0.17 

± 0.16 

0.48 

± 

0.21 

± 

0.35 

0.88 

± 0.12 

± 0.24 

0.070 - 0.100 

0.320 

± 0.004 

± 0.013 

0.204 ± 0.076 ± 0.072 

0.13 

± 

0.20 

± 

0.33 

0.03 

± 0.16 

± 0.51 


Table 4: Distributions for the thrust major (Tmaj.), thrust minor (Tmin.) and aplanarity (A) measured by OPAL at centre-of-mass 
energies ^/s = 91, 130-136, 161-183, and 189-209 GeV. The hrst uncertainty is statistical, while the second is systematic. 




s 

R ! { S ) at 91 GeV 

R ' ( S ) at 133 GeV 

R ' ( S ) at 177 GeV 

R ' ( S ) at 197 GeV 

0.00 - 0.02 

18.06 

± 0.04 

± 0.42 

23.63 

± 0.95 

± 1.81 

25.66 

± 0.70 

± 0.56 

25.67 ± 0.32 ± 0.37 

0.02 - 0.04 

10.57 

± 0.03 

± 0.15 

8.15 

± 0.72 

± 0.77 

6.86 

± 0.46 

± 0.29 

7.36 ± 0.20 ± 0.30 

0.04 - 0.06 

5.183 

± 0.024 

± 0.074 

4.07 

± 0.53 

± 0.28 

3.83 

± 0.35 

± 0.49 

3.65 ± 0.15 ± 0.18 

0.06 - 0.12 

2.409 

± 0.009 

± 0.058 

1.65 

± 0.22 

± 0.13 

2.21 

± 0.15 

± 0.25 

1.99 ± 0.06 ± 0.13 

0.12 - 0.20 

0.988 

± 0.005 

± 0.033 

0.93 

± 0.12 

± 0.22 

0.74 

± 0.09 

± 0.12 

0.784 ± 0.041 ± 0.070 

0.20 - 0.30 

0.466 

± 0.003 

± 0.014 

0.377 ± 0.068 ± 0.048 

0.338 ± 0.070 ± 0.073 

0.379 ± 0.033 ± 0.044 

0.30 - 0.50 

0.2001 ± 0.0015 ± 0.0100 

0.293 ± 0.034 ± 0.033 

0.138 ± 0.040 ± 0.028 

0.175 ± 0.023 ± 0.047 

0.50 - 0.70 

0.0614 ± 0.0008 ± 0.0026 

0.068 ± 0.016 ± 0.013 

0.065 ± 0.030 ± 0.036 

0.044 ± 0.019 ± 0.028 

0 

H ( O ) at 91 GeV 

R ' ( 0 ) at 133 GeV 

R ' ( 0 ) at 177 GeV 

R ' ( 0 ) at 197 GeV 

0.00 - 0.05 

9.934 

± 0.016 

± 0.142 

9.71 

± 0.38 

± 0.20 

10.07 

± 0.26 

± 0.32 

9.976 ± 0.111 ± 0.212 

0.05 - 0.10 

4.576 

± 0.013 

± 0.039 

4.52 

± 0.31 

± 0.24 

4.33 

± 0.21 

± 0.36 

4.188 ± 0.089 ± 0.047 

0.10 - 0.15 

2.276 

± 0.010 

± 0.041 

2.17 

± 0.25 

± 0.13 

2.12 

± 0.18 

± 0.30 

2.255 ± 0.073 ± 0.052 

0.15 - 0.20 

1.307 

± 0.008 

± 0.030 

1.20 

± 0.19 

± 0.13 

1.34 

± 0.14 

± 0.17 

1.287 ± 0.062 ± 0.054 

0 . 20 - 0.25 

0.810 

± 0.006 

± 0.017 

0.82 

± 0.15 

± 0.07 

0.87 

± 0.12 

± 0.07 

0.837 ± 0.053 ± 0.064 

0.25 - 0.30 

0.490 

± 0.005 

± 0.024 

0.65 

± 0.11 

± 0.13 

0.47 

± 0.10 

± 0.16 

0.596 ± 0.046 ± 0.059 

0.30 - 0.40 

0.2361 ± 0.0024 ± 0.0071 

0.318 ± 0.062 ± 0.045 

0.32 

± 0.06 

± 0.10 

0.307 ± 0.026 ± 0.052 

0.40 - 0.50 

0.0521 ± 0.0011 ± 0.0028 

0.150 ± 0.028 ± 0.059 

0.073 ± 0.030 ± 0.043 

0.108 ± 0.016 ± 0.057 

Ml 

R ' (Ml) at 91 GeV 

R ' (Ml) at 133 GeV 

R ' (Ml) at 177 GeV 

R ' (Ml) at 197 GeV 

0.00 - 0.04 

0.124 

± 0.002 

± 0.038 

0.292 ± 0.067 ± 0.080 

0.312 ± 0.056 ± 0.050 

0.430 ± 0.027 ± 0.033 

0.04 - 0.06 

1.112 

± 0.008 

± 0.096 

2.29 

± 0.29 

± 0.57 

4.81 

± 0.29 

± 0.39 

4.51 ± 0.14 ± 0.36 

0.06 - 0.08 

3.72 

± 0.02 

± 0.56 

7.5 

± 0.6 

± 1.4 

9.98 

± 0.51 

± 0.44 

11.41 ± 0.24 ± 0.31 

0.08 - 0.10 

8.22 

± 0.03 

± 0.57 

11.1 

± 0.8 

± 1.7 

11.89 

± 0.59 

± 1.02 

11.72 ± 0.26 ± 0.63 

0.10 - 0.12 

10.73 

± 0.03 

± 0.36 

10.3 

± 0.8 

± 1.1 

8.47 

± 0.54 

± 1.30 

7.95 ± 0.22 ± 0.26 

0.12 - 0.14 

9.11 

± 0.03 

± 0.55 

6.2 

± 0.7 

± 1.0 

5.17 

± 0.43 

± 0.49 

4.68 ± 0.17 ± 0.24 

0.14 - 0.16 

6.23 

± 0.03 

± 0.23 

4.7 

± 0.6 

± 1.1 

2.99 

± 0.36 

± 0.29 

3.08 ± 0.15 ± 0.12 

0.16 - 0.20 

3.150 

± 0.015 

± 0.086 

1.85 

± 0.28 

± 0.27 

1.83 

± 0.20 

± 0.16 

1.69 ± 0.08 ± 0.12 

0.20 - 0.24 

1.343 

± 0.009 

± 0.031 

1.03 

± 0.18 

± 0.15 

0.53 

± 0.15 

± 0.29 

0.779 ± 0.067 ± 0.059 

0.24 - 0.30 

0.450 

± 0.004 

± 0.016 

0.461 ± 0.085 ± 0.089 

0.36 

± 0.09 

± 0.20 

0.296 ± 0.052 ± 0.057 

0.30 - 0.40 

0.0609 ± 0.0013 ± 0.0041 

0.038 ± 0.025 ± 0.024 

0.05 

± 0.04 

± 0.13 

- 0.008 ± 0.061 ± 0.086 


Table 5: Distributions for the sphericity (S'), oblateness (O) and light jet mass (Ml) measured by OPAL at centre-of-mass 
energies ^/s = 91, 130-136, 161-183, and 189-209 GeV. The hrst uncertainty is statistical, while the second is systematic. 




Sn 

B ! ( Bn ) at 91 GeV 

R ' (- Bn ) at 133 GeV 

R ' (- Bn ) at 177 GeV 

R ' (- Bn ) at 197 GeV 

0.000 - 

0.010 

0.698 

± 0.020 

± 0.081 

1.83 

± 0.71 

± 1.00 

4.86 

± 0.69 

± 0.81 

6.33 

± 0.34 

± 0.29 

0 . 010 - 

0.015 

7.96 

± 0.07 

± 0.20 

30.0 

± 2.8 

± 1.6 

38.2 

± 2.2 

± 2.4 

42.4 

± 1.1 

± 1.6 

0 . 015 - 

0.020 

21.57 

± 0.11 

± 0.36 

37.7 

± 3.2 

± 4.9 

50.9 

± 2.3 

± 2.3 

45.1 

± 1.1 

± 1.7 

0.020 - 

0.025 

32.46 

± 0.12 

± 0.47 

34.1 

± 3.2 

± 6.9 

29.8 

± 1.8 

± 3.0 

27.8 

± 0.8 

± 1.6 

0.025 - 

0.030 

33.65 

± 0.12 

± 0.73 

26.3 

± 2.5 

± 3.7 

19.5 

± 1.4 

± 1.5 

16.6 

± 0.6 

± 0.8 

0.030 - 

0.035 

28.42 

± 0.11 

± 0.17 

15.6 

± 2.0 

± 1.0 

11.4 

± 1.2 

± 1.3 

13.0 

± 0.5 

± 1.3 

0.035 - 

0.040 

19.94 

± 0.09 

± 0.31 

10.3 

± 1.7 

± 1.2 

8.1 

± 1.0 

± 1.9 

8.45 

± 0.41 

± 0.66 

0.040 - 

0.050 

11.47 

± 0.05 

± 0.41 

8.2 

± 0.9 

± 2.4 

6.81 

± 0.59 

± 0.38 

5.53 

± 0.26 

± 0.40 

0.050 - 

0.060 

6.053 

± 0.034 

± 0.084 

4.15 

± 0.68 

± 0.85 

3.22 

± 0.49 

± 0.70 

3.80 

± 0.22 

± 0.36 

0.060 - 

0.080 

2.914 

± 0.017 

± 0.046 

2.81 

± 0.35 

± 0.23 

1.85 

± 0.31 

± 0.53 

1.58 

± 0.14 

± 0.16 

0.080 - 

0.120 

0.835 

± 0.006 

± 0.017 

0.68 

± 0.14 

± 0.12 

0.52 

± 0.18 

± 0.42 

0.73 

± 0.10 

± 0.21 

0 . 120 - 

0.170 

0.1211 ± 0.0021 ± 0.0043 

0.094 

± 0.042 

± 0.035 

0.07 

± 0.13 

± 0.27 

0.33 

± 0.23 

± 0.37 

D 

R ' ( D ) at 91 GeV 

R ' { D ) at 133 GeV 

R ' ( D ) at 177 GeV 

R ' ( D ) at 197 GeV 

0.001 - 

0.005 

23.07 

± 0.12 

± 1.16 

43.9 

± 3.9 

± 7.9 

63.3 

± 2.9 

± 3.0 

61.6 

± 1.3 

± 1.3 

0.005 - 

0.010 

29.68 

± 0.12 

± 0.34 

40.9 

± 2.9 

± 3.0 

29.5 

± 1.8 

± 2.5 

28.4 

± 0.8 

± 1.4 

0 . 010 - 

0.015 

21.76 

± 0.10 

± 0.25 

20.2 

± 2.2 

± 1.3 

17.4 

± 1.4 

± 1.8 

16.71 

± 0.59 

± 0.50 

0.015 - 

0.020 

15.62 

± 0.09 

± 0.40 

10.5 

± 1.9 

± 2.4 

9.0 

± 1.2 

± 1.8 

12.18 

± 0.50 

± 0.27 

0.020 - 

0.030 

10.85 

± 0.05 

± 0.12 

8.4 

± 1.1 

± 1.1 

8.07 

± 0.71 

± 0.60 

7.99 

± 0.29 

± 0.20 

0.030 - 

0.045 

6.651 

± 0.032 

± 0.097 

4.37 

± 0.72 

± 0.93 

5.69 

± 0.46 

± 0.61 

5.32 

± 0.19 

± 0.23 

0.045 - 

0.070 

4.031 

± 0.019 

± 0.031 

3.31 

± 0.43 

± 0.26 

2.81 

± 0.28 

± 0.22 

3.02 

± 0.12 

± 0.11 

0.070 - 

0.100 

2.400 

± 0.014 

± 0.035 

2.24 

± 0.30 

± 0.32 

1.35 

± 0.20 

± 0.22 

1.76 

± 0.09 

± 0.12 

0 . 100 - 

0.150 

1.424 

± 0.008 

± 0.027 

1.06 

± 0.17 

± 0.19 

1.27 

± 0.13 

± 0.19 

0.928 ± 0.055 ± 0.075 

0.150 - 

0.250 

0.672 

± 0.004 

± 0.024 

0.66 

± 0.08 

± 0.12 

0.36 

± 0.07 

± 0.11 

0.446 ± 0.032 ± 0.060 

0.250 - 

0.500 

0.1660 ± 0.0012 ± 0.0020 

0.122 

± 0.022 

± 0.024 

0.133 ± 0.032 ± 0.051 

0.116 ± 0.019 ± 0.043 

0.500 - 

1.000 

0.0141 ± 0.0002 ± 0.0006 

0.0042 ± 0.0036 ± 0.0020 

0.033 ± 0.017 ± 0.027 

0.017 ± 0.011 ± 0.021 


Table 6: Distributions for the narrow jet broadening (-Bn) and .D-parameter measured by OPAL at centre-of-mass energies 
y/s = 91, 130-136, 161-183, and 189-209 GeV. The first uncertainty is statistical, while the second is systematic. 




n 

((1 - Tf ) at 91 GeV 

(( 1 - r )") at 133 GeV 

(( 1 - r )") at 177 GeV 

((1 - r )") at 197 GeV 

1 

( 6.671 

± 0.017 

± 0.066 ) 

10-2 

( 6.20 

± 0.26 

± 0.32 ) 

10-2 

( 5.60 

± 0.17 

± 0.22 ) 

10-2 

( 5.60 

± 0.08 

± 0.15 ) 

10-2 

2 

( 8.12 

± 0.04 

± 0.13 ) 

10-2 

( 7.84 

± 0.64 

± 0.53 ) 

10-2 

( 6.78 

± 0.48 

± 0.72 ) 

10-2 

( 6.76 

± 0.22 

± 0.41 ) 

10-5 

3 

( 1.457 

± 0.012 

± 0.025 ) 

10-2 

( 1.43 

± 0.17 

± 0.11 ) 

10-2 

( 1.26 

± 0.15 

± 0.24 ) 

10-2 

( 1.22 

± 0.07 

± 0.12 ) 

10-5 

4 

( 3.236 

± 0.037 

± 0.053 ) 

10-2 

( 3.13 

± 0.48 

± 0.33 ) 

10-2 

( 2.96 

± 0.47 

± 0.87 ) 

10-2 

( 2.71 

± 0.25 

± 0.41 ) 

10-5 

5 

( 8.13 

± 0.12 

± 0.13 ) 

10-2 

( 7.6 

± 1.5 

± 1.1 ) 

10-2 

( 8.0 

± 1.5 

± 3.3 ) 

10-5 

( 6.7 

± 0.9 

± 1.6 ) 

10-5 

n 

( MS ) at 91 GeV 

( Mg ) at 133 GeV 

( Mg ) at 177 GeV 

( Mg ) at 197 GeV 

1 

( 2.106 

± 0.003 

± 0.018 ) 

10-2 

( 1.992 ± 0.043 ± 0 . 063 ) 

10-2 

( 1.887 ± 0.028 ± 0 . 038 ) 

10-2 

( 1.877 ± 0.013 ± 0 . 030 ) 

10-5 

2 

( 5.235 

± 0.014 

± 0.086 ) 

10-2 

( 4.96 

± 0.22 

± 0.27 ) 

10-2 

( 4.50 

± 0.15 

± 0.22 ) 

10-2 

( 4.51 

± 0.07 

± 0.17 ) 

10-2 

3 

( 1.523 

± 0.006 

± 0.036 ) 

10-2 

( 1.50 

± 0.11 

± 0.10 ) 

10-2 

( 1.31 

± 0.07 

± 0.11 ) 

10-2 

( 1.339 ± 0.034 ± 0 . 086 ) 

10-2 

4 

( 5.06 

± 0.03 

± 0.15 ) 

10-2 

( 5.26 

± 0.50 

± 0.44 ) 

10-2 

( 4.39 

± 0.33 

± 0.53 ) 

10-2 

( 4.63 

± 0.17 

± 0.43 ) 

10-5 

5 

( 1.862 

± 0.015 

± 0.063 ) 

10-2 

( 2.04 

± 0.24 

± 0.21 ) 

10-2 

( 1.62 

± 0.16 

± 0.25 ) 

10-2 

( 1.78 

± 0.09 

± 0.22 ) 

10-5 

n 

( C ") at 91 GeV 

( G ") at 133 GeV 

( C ") at 177 GeV 

( C ") at 197 GeV 

1 

( 2.646 

± 0.005 

± 0.022 ) 

10-2 

( 2.42 

± 0.08 

± 0.11 ) 

10-2 

( 2.210 ± 0.053 ± 0 . 064 ) 

10-1 

( 2.214 ± 0.024 ± 0 . 042 ) 

10-5 

2 

( 1.049 

± 0.004 

± 0.015 ) 

10-2 

( 9.70 

± 0.59 

± 0.64 ) 

10-2 

( 8.45 

± 0.41 

± 0.54 ) 

10-2 

( 8.52 

± 0.19 

± 0.31 ) 

10-2 

3 

( 5.412 

± 0.030 

± 0.094 ) 

10-2 

( 5.13 

± 0.42 

± 0.37 ) 

10-2 

( 4.38 

± 0.31 

± 0.44 ) 

10-2 

( 4.39 

± 0.15 

± 0.22 ) 

10-2 

4 

( 3.244 

± 0.023 

± 0.060 ) 

10-2 

( 3.10 

± 0.31 

± 0.23 ) 

10-2 

( 2.65 

± 0.25 

± 0.37 ) 

10-2 

( 2.62 

± 0.13 

± 0.16 ) 

10-2 

5 

( 2.125 

± 0.018 

± 0.040 ) 

10-2 

( 2.02 

± 0.23 

± 0.15 ) 

10-2 

( 1.76 

± 0.20 

± 0.33 ) 

10-2 

( 1.69 

± 0.11 

± 0.15 ) 

10-2 

n 

( B ?) at 91 GeV 

{ BJ }.) at 133 GeV 

( BJ ^) at 177 GeV 

( B ?) at 197 GeV 

1 

( 1.0909 ± 0.0016 ± 0 . 0068 ) 

10-2 

( 1.007 ± 0.026 ± 0 . 036 ) 

10-2 

( 9.28 

± 0.17 

± 0.23 ) 

10-2 

( 9.25 

± 0.08 

± 0.15 ) 

10-2 

2 

( 1.548 

± 0.005 

± 0.018 ) 

10-2 

( 1.428 ± 0.075 ± 0 . 083 ) 

10-2 

( 1.244 ± 0.050 ± 0 . 071 ) 

10-2 

( 1.245 ± 0.023 ± 0 . 045 ) 

10-2 

3 

( 2.738 

± 0.014 

± 0.042 ) 

10-2 

( 2.60 

± 0.20 

± 0.18 ) 

10-2 

( 2.19 

± 0.14 

± 0.22 ) 

10-2 

( 2.19 

± 0.07 

± 0.12 ) 

10-5 

4 

( 5.684 

± 0.039 

± 0.096 ) 

10-2 

( 5.53 

± 0.56 

± 0.42 ) 

10-2 

( 4.62 

± 0.44 

± 0.71 ) 

10-2 

( 4.57 

± 0.22 

± 0.34 ) 

10-5 

5 

( 1.317 

± 0.012 

± 0.023 ) 

10-2 

( 1.29 

± 0.16 

± 0.11 ) 

10-2 

( 1.10 

± 0.14 

± 0.25 ) 

10-2 

( 1.06 

± 0.07 

± 0.11 ) 

10-5 

n 

( B ^) at 91 GeV 

( B ^) at 133 GeV 

( B ^) at 177 GeV 

( B ^) at 197 GeV 

1 

( 7.361 

± 0.013 

± 0.065 ) 

10-2 

( 7.04 

± 0.21 

± 0.28 ) 

10-2 

( 6.63 

± 0.14 

± 0.21 ) 

10-2 

( 6.65 

± 0.06 

± 0.15 ) 

10-2 

2 

( 7.70 

± 0.03 

± 0.13 ) 

10-2 

( 7.71 

± 0.48 

± 0.50 ) 

10-2 

( 6.94 

± 0.33 

± 0.50 ) 

10-2 

( 7.09 

± 0.15 

± 0.36 ) 

10-5 

3 

( 1.055 

± 0.006 

± 0.024 ) 

10-2 

( 1.14 

± 0.10 

± 0.09 ) 

10-2 

( 9.8 

± 0.7 

± 1.1 ) 

10-2 

( 1.022 ± 0.037 ± 0 . 084 ) 

10-5 

4 

( 1.723 

± 0.014 

± 0.045 ) 

10-2 

( 2.01 

± 0.24 

± 0.20 ) 

10-2 

( 1.63 

± 0.17 

± 0.26 ) 

10-2 

( 1.76 

± 0.09 

± 0.21 ) 

10-5 

5 

( 3.150 

± 0.032 

± 0.088 ) 

10-2 

( 3.93 

± 0.58 

± 0.47 ) 

10-2 

( 3.00 

± 0.41 

± 0.61 ) 

10-5 

( 3.36 

± 0.23 

± 0.52 ) 

10-5 

77 

( 2 / 2 °") at 91 GeV 

{ yfs ") at 133 GeV 

{ yg ") at 177 GeV 

( j / D ”) at 197 GeV 

1 

( 2.059 

± 0.011 

± 0.051 ) 

10-2 

( 2.14 

± 0.16 

± 0.12 ) 

10-2 

( 1.87 

± 0.12 

± 0.19 ) 

10-2 

( 1.91 

± 0.05 

± 0.13 ) 

10-2 

2 

( 1.951 

± 0.020 

± 0.047 ) 

10-2 

( 2.21 

± 0.31 

± 0.18 ) 

10-2 

( 1.74 

± 0.23 

± 0.42 ) 

10-2 

( 1.86 

± 0.12 

± 0.27 ) 

10-5 

3 

( 3.053 

± 0.046 

± 0.064 ) 

10-2 

( 3.60 

± 0.76 

± 0.54 ) 

10-2 

( 2.7 

± 0.5 

± 1.0 ) 

10-2 

( 2.90 

± 0.31 

± 0.61 ) 

10-5 

4 

( 5.89 

± 0.12 

± 0.15 ) 

10-2 

( 7.2 

± 2.0 

± 1.7 ) 

10-2 

( 5.1 

± 1.0 

± 2.8 ) 

10-5 

( 5.5 

± 0.8 

± 1.5 ) 

10-5 

5 

( 1.270 

± 0.031 

± 0.043 ) 

10-2 

( 1.61 

± 0.53 

± 0.54 ) 

10-2 

( 1.14 

± 0.20 

± 0.80 ) 

10-5 

( 1.12 

± 0.21 

± 0.37 ) 

10-5 


Table 7: Moments of the (1 — T), Mh, C, B-r, i?w and distributions measured by OPAL at 91, 130-136, 161-183, and 
189-209 GeV. The first uncertainty is statistical, while the second is systematic. 





























n 

( r , Saj .) at 91 GeV 

at 133 GeV 

( r , Saj .) at 177 GeV 

( r , Saj .) at 197 GeV 

1 

( 1.757 

± 0.003 ± 0.013 ) 

10-1 

( 1.671 ± 0.048 ± 0 . 059 ) 

10-1 

( 1.546 ± 0.031 ± 0 . 042 ) 

10-1 

( 1.550 ± 0.014 ± 0 . 031 ) 

10-5 

2 

( 4.244 

± 0.015 ± 0.059 ) 

10-2 

( 4.20 

± 0.24 

± 0.24 ) 

10-2 

( 3.65 

± 0.16 

± 0.22 ) 

10-2 

( 3.71 

± 0.08 ± 0.17 ) 

10-5 

3 

( 1.321 

± 0.007 ± 0.024 ) 

10-2 

( 1.39 

± 0.12 

± 0.10 ) 

10-2 

( 1.15 

± 0.08 

± 0.12 ) 

10-2 

( 1.190 ± 0.042 ± 0 . 082 ) 

10-5 

4 

( 4.88 

± 0.04 ± 0.10 ) 

10-3 

( 5.43 

± 0.59 

± 0.42 ) 

10-2 

( 4.33 

± 0.45 

± 0.63 ) 

10-5 

( 4.55 

± 0.24 ± 0.42 ) 

10-5 

5 

( 2.019 

± 0.019 ± 0.044 ) 

10-3 

( 2.33 

± 0.30 

± 0.20 ) 

10-2 

( 1.81 

± 0.25 

± 0.35 ) 

10-5 

( 1.94 

± 0.14 ± 0.23 ) 

10-5 

n 

at 91 GeV 

( T ";,,) at 133 GeV 

at 177 GeV 

( r , Si ,,) at 197 GeV 

1 

( 9.553 

± 0.012 ± 0.030 ) 

10-2 

( 7.98 

± 0.17 

± 0.33 ) 

10-2 

( 7.25 

± 0.11 

± 0.12 ) 

10-2 

( 7.106 ± 0.049 ± 0 . 080 ) 

10-5 

2 

( 1.1160 ± 0.0034 ± 0 . 0061 ) 

10-2 

( 8.15 

± 0.42 

± 0.63 ) 

10-2 

( 7.21 

± 0.30 

± 0.38 ) 

10-5 

( 6.97 

± 0.13 ± 0.26 ) 

10-5 

3 

( 1.614 

± 0.010 ± 0.012 ) 

10-3 

( 1.06 

± 0.10 

± 0.11 ) 

10-2 

( 1.04 

± 0.10 

± 0.16 ) 

10-5 

( 9.61 

± 0.41 ± 0.90 ) 

10-5 

4 

( 2.879 

± 0.031 ± 0.032 ) 

10-3 

( 1.72 

± 0.27 

± 0.24 ) 

10-2 

( 2.26 

± 0.39 

± 0.87 ) 

10-5 

( 1.80 

± 0.16 ± 0.38 ) 

10-5 

5 

( 6.20 

± 0.11 ± 0.10 ) 

10-3 

( 3.30 

± 0.74 

± 0.64 ) 

10 -® 

( 7.5 

± 1.6 

± 5.5 ) 

10-5 

( 4.2 

± 0.7 ± 2.0 ) 

10-5 

n 

( S ") at 91 GeV 

( 5 "*) at 133 GeV 

( S ") at 177 GeV 

( S ") at 197 GeV 

1 

( 7.42 

± 0.03 ± 0.20 ) 

10-2 

( 7.31 

± 0.48 

± 0.45 ) 

10-2 

( 6.10 

± 0.33 

± 0.46 ) 

10-2 

( 6.00 

± 0.15 ± 0.30 ) 

10-5 

2 

( 1.706 

± 0.014 ± 0.069 ) 

10-2 

( 1.89 

± 0.22 

± 0.14 ) 

10-2 

( 1.52 

± 0.21 

± 0.30 ) 

10-2 

( 1.41 

± 0.10 ± 0.10 ) 

10-5 

3 

( 6.43 

± 0.08 ± 0.30 ) 

10-3 

( 7.4 

± 1.2 

± 0.8 ) 

10-2 

( 6.8 

± 1.4 

± 2.4 ) 

10-5 

( 5.42 

± 0.77 ± 0.74 ) 

10-5 

4 

( 3.01 

± 0.05 ± 0.15 ) 

10-3 

( 3.43 

± 0.69 

± 0.51 ) 

10-2 

( 3.9 

± 0.9 

± 1.9 ) 

10-5 

( 2.56 

± 0.59 ± 0.73 ) 

10-5 

5 

( 1.586 

± 0.033 ± 0.088 ) 

10-3 

( 1.77 

± 0.42 

± 0.35 ) 

10-2 

( 2.5 

± 0.6 

± 1.5 ) 

10-5 

( 1.38 

± 0.44 ± 0.64 ) 

10-5 

n 

( 0 ") at 91 GeV 

( O ") at 133 GeV 

( 0 ") at 177 GeV 

( 0 ") at 197 GeV 

1 

( 8.03 

± 0.02 ± 0.10 ) 

10-2 

( 8.74 

± 0.38 

± 0.29 ) 

10-2 

( 8.21 

± 0.26 

± 0.32 ) 

10-2 

( 8.39 

± 0.12 ± 0.31 ) 

10-5 

2 

( 1.295 

± 0.007 ± 0.030 ) 

10-2 

( 1.62 

± 0.13 

± 0.10 ) 

10-2 

( 1.40 

± 0.09 

± 0.11 ) 

10-2 

( 1.52 

± 0.05 ± 0.12 ) 

10-5 

3 

( 3.031 

± 0.026 ± 0.086 ) 

10-3 

( 4.28 

± 0.50 

± 0.40 ) 

10-2 

( 3.39 

± 0.37 

± 0.41 ) 

10-5 

( 4.00 

± 0.20 ± 0.54 ) 

10-5 

4 

( 8.67 

± 0.10 ± 0.27 ) 

10-3 

( 1.34 

± 0.20 

± 0.18 ) 

10-2 

( 9.9 

± 1.5 

± 1.7 ) 

10-5 

( 1.29 

± 0.09 ± 0.24 ) 

10-5 

5 

( 2.803 

± 0.041 ± 0.096 ) 

10-3 

( 4.64 

± 0.81 

± 0.80 ) 

10-2 

( 3.20 

± 0.65 

± 0.76 ) 

10-5 

( 4.6 

± 0.4 ± 1.1 ) 

10-5 

n 

(MO at 91 GeV 

( MO ) at 133 GeV 

{ M £) at 177 GeV 

{ Ml ) at 197 GeV 

1 

( 1.3111 ± 0.0014 ± 0 . 0050 ) 

10-3 

( 1.162 ± 0.022 ± 0 . 026 ) 

10-1 

( 1.059 ± 0.014 ± 0 . 013 ) 

10-1 

( 1.047 ± 0.006 ± 0 . 011 ) 

10-1 

2 

( 1.945 

± 0.005 ± 0.016 ) 

10-2 

( 1.590 ± 0.068 ± 0 . 054 ) 

10-2 

( 1.337 ± 0.041 ± 0 . 044 ) 

10-2 

( 1.311 ± 0.018 ± 0 . 032 ) 

10-5 

3 

( 3.265 

± 0.013 ± 0.038 ) 

10-3 

( 2.56 

± 0.18 

± 0.08 ) 

10-2 

( 2.05 

± 0.12 

± 0.15 ) 

10-5 

( 1.970 ± 0.049 ± 0 . 091 ) 

10-5 

4 

( 6.176 

± 0.039 ± 0.086 ) 

10-3 

( 4.75 

± 0.49 

± 0.13 ) 

10-2 

( 3.81 

± 0.39 

± 0.54 ) 

10-5 

( 3.51 

± 0.15 ± 0.27 ) 

10-5 

5 

( 1.304 

± 0.012 ± 0.020 ) 

10-3 

( 10.0 

± 1.4 

± 0.5 ) 

10 -® 

( 8.6 

± 1.4 

± 2.2 ) 

10-5 

( 7.23 

± 0.47 ± 0.85 ) 

10-5 

n 

( B ^) at 91 GeV 

( BJ )) at 133 GeV 

( S ( 5 ) at 177 GeV 

( Bl ) at 197 GeV 

1 

( 3.5456 ± 0.0054 ± 0 . 0089 ) 

10-2 

( 3.029 ± 0.078 ± 0 . 077 ) 

10-2 

( 2.649 ± 0.048 ± 0 . 043 ) 

10-2 

( 2.600 ± 0.022 ± 0 . 029 ) 

10-5 

2 

( 1.6369 ± 0.0062 ± 0 . 0058 ) 

10-3 

( 1.302 ± 0.081 ± 0 . 054 ) 

10-2 

( 1.062 ± 0.055 ± 0 . 062 ) 

10-5 

( 1.042 ± 0.025 ± 0 . 040 ) 

10-5 

3 

( 1.0016 ± 0.0073 ± 0 . 0042 ) 

10-3 

( 7.80 

± 0.90 

± 0.39 ) 

10-5 

( 6.7 

± 0.8 

± 1.1 ) 

10-5 

( 6.41 

± 0.37 ± 0.59 ) 

10-5 

4 

( 7.900 

± 0.094 ± 0.072 ) 

10 -® 

( 6.0 

± 1.1 

± 0.4 ) 

10-5 

( 6.3 

± 1.4 

± 2.6 ) 

10-5 

( 5.5 

± 0.7 ± 1.0 ) 

10-5 

5 

( 7.57 

± 0.13 ± 0.14 ) 

10-2 

( 5.6 

± 1.4 

± 0.5 ) 

10 -^ 

( 8.6 

± 2.2 

± 6.2 ) 

10 -^ 

( 5.8 

± 1.4 ± 2.0 ) 

10 -^ 


Table 8: Moments of the Tmaj., T^in., S, O, Mi^ and -Bn distributions measured by OPAL at 91, 130-136, 161-183, and 189 
209 GeV. The first uncertainty is statistical, while the second is systematic. 





























T 

Mu 


Sw 

C 

2/23 

Weighted 

mean 

as (91 GeV) 

0.1231 

0.1192 

0.1224 

0.1146 

0.1178 

0.1197 

0.1192 

Statistical error 

±0.0002 

±0.0002 

±0.0002 

±0.0002 

±0.0002 

±0.0002 

±0.0002 

Experimental syst. 

±0.0011 

±0.0004 

±0.0007 

±0.0011 

±0.0007 

±0.0010 

±0.0008 

HERWIG hadr. corr. 

-0.0018 

±0.0013 

-0.0023 

-0.0014* 

-0.0031* 

-0.0024* 

-0.0015* 

ARIADNE hadr. corr. 

±0.0031* 

±0.0021* 

±0.0024* 

±0.0007 

±0.0029 

-0.0008 

±0.0013 

Hadronization error 

±0.0031 

±0.0021 

±0.0024 

±0.0014 

±0.0031 

±0.0024 

±0.0015 

Theory error 

±0.0054 

±0.0043 

±0.0064 

±0.0053 

±0.0055 

±0.0033 

±0.0047 

Weight 

0.12 

0.21 

0.11 

0.16 

0.12 

0.28 

- 

as(133 GeV) 

0.1158 

0.1088 

0.1102 

0.1051 

0.1056 

0.1109 

0.1092 

Statistical error 

±0.0041 

±0.0038 

±0.0039 

±0.0032 

±0.0043 

±0.0031 

±0.0032 

Experimental syst. 

±0.0054 

±0.0025 

±0.0063 

±0.0047 

±0.0052 

±0.0047 

±0.0044 

HERWIG hadr. corr. 

-0.0005 

±0.0020* 

-0.0016 

-0.0005 

-0.0014 

-0.0012* 

-0.0002 

ARIADNE hadr. corr. 

±0.0024* 

±0.0018 

±0.0020* 

±0.0007* 

±0.0026* 

-0.0005 

±0.0012* 

Hadronization error 

±0.0024 

±0.0020 

±0.0020 

±0.0007 

±0.0026 

±0.0012 

±0.0012 

Theory error 

±0.0047 

±0.0037 

±0.0055 

±0.0046 

±0.0047 

±0.0028 

±0.0040 

Weight 

0.12 

0.25 

0.08 

0.17 

0.12 

0.25 

- 

Os (177 GeV) 

0.1105 

0.1068 

0.1085 

0.1016 

0.1062 

0.1070 

0.1063 

Statistical error 

±0.0030 

±0.0027 

±0.0026 

±0.0023 

±0.0029 

±0.0022 

±0.0021 

Experimental syst. 

±0.0028 

±0.0033 

±0.0033 

±0.0021 

±0.0024 

±0.0031 

±0.0026 

HERWIG hadr. corr. 

-0.0001 

±0.0021* 

-0.0010 

-0.0002 

-0.0006 

-0.0005* 

-0.0001 

ARIADNE hadr. corr. 

±0.0021* 

±0.0015 

±0.0016* 

±0.0006* 

±0.0020* 

-0.0003 

±0.0010* 

Hadronization error 

±0.0021 

±0.0021 

±0.0016 

±0.0006 

±0.0020 

±0.0005 

±0.0010 

Theory error 

±0.0042 

±0.0033 

±0.0050 

±0.0041 

±0.0042 

±0.0024 

±0.0036 

Weight 

0.12 

0.14 

0.10 

0.20 

0.14 

0.30 

- 

as(197 GeV) 

0.1130 

0.1063 

0.1118 

0.1031 

0.1077 

0.1071 

0.1075 

Statistical error 

±0.0013 

±0.0011 

±0.0018 

±0.0009 

±0.0012 

±0.0009 

±0.0010 

Experimental syst. 

±0.0013 

±0.0016 

±0.0011 

±0.0010 

±0.0012 

±0.0010 

±0.0009 

HERWIG hadr. corr. 

-0.0001 

±0.0021* 

-0.0007 

-0.0001 

-0.0004 

-0.0003* 

0.0000 

ARIADNE hadr. corr. 

±0.0019* 

±0.0014 

±0.0013* 

±0.0006* 

±0.0018* 

-0.0003 

±0.0007* 

Hadronization error 

±0.0019 

±0.0021 

±0.0013 

±0.0006 

±0.0018 

±0.0003 

±0.0007 

Theory error 

±0.0040 

±0.0032 

±0.0048 

±0.0040 

±0.0041 

±0.0023 

±0.0033 

Weight 

0.11 

0.14 

0.08 

0.16 

0.11 

0.40 

- 


Table 9: Measurements of ag using event shape distributions in four ranges of c.m. en¬ 
ergy: at 91 GeV, 133 GeV, 161-183 GeV (denoted 177 GeV) and 189-209 GeV (denoted 
197 GeV). The hadronization error is taken to be the larger of the effects observed using 
HERWIG and ARIADNE; in each case this is denoted by an asterisk. The weights and 
weighted mean are described in the text. 
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T 

Mu 

Bt 

Bw 

c 

y?3 

Weighted 

mean 

as(161 GeV) 

0.1103 

0.1064 

0.1051 

0.1010 

0.1042 

0.1042 

0.1046 

Statistical error 

±0.0069 

±0.0063 

±0.0062 

±0.0053 

±0.0068 

±0.0051 

±0.0051 

Experimental syst. 

±0.0042 

±0.0042 

±0.0045 

±0.0035 

±0.0056 

±0.0033 

±0.0034 

HERWIG hadr. corr. 

-0.0002 

±0.0021* 

-0.0013 

-0.0002 

-0.0009 

-0.0009* 

-0.0002 

ARIADNE hadr. corr. 

±0.0022* 

±0.0015 

±0.0018* 

±0.0007* 

±0.0022* 

-0.0005 

±0.0009* 

Hadronization error 

±0.0022 

±0.0021 

±0.0018 

±0.0007 

±0.0022 

±0.0009 

±0.0009 

Theory error 

±0.0043 

±0.0034 

±0.0052 

±0.0043 

±0.0044 

±0.0025 

±0.0036 

Weight 

0.12 

0.14 

0.06 

0.21 

0.09 

0.39 

- 

as(172 GeV) 

0.1095 

0.1043 

0.1016 

0.0950 

0.1039 

0.1038 

0.1011 

Statistical error 

±0.0077 

±0.0069 

±0.0070 

±0.0060 

±0.0075 

±0.0056 

±0.0055 

Experimental syst. 

±0.0076 

±0.0078 

±0.0037 

±0.0029 

±0.0064 

±0.0046 

±0.0039 

HERWIG hadr. corr. 

-0.0002 

±0.0022* 

-0.0012 

-0.0002 

-0.0007 

-0.0007* 

-0.0003 

ARIADNE hadr. corr. 

±0.0023* 

±0.0017 

±0.0018* 

±0.0008* 

±0.0022* 

-0.0003 

±0.0008* 

Hadronization error 

±0.0023 

±0.0022 

±0.0018 

±0.0008 

±0.0022 

±0.0007 

±0.0008 

Theory error 

±0.0042 

±0.0033 

±0.0050 

±0.0041 

±0.0042 

±0.0024 

±0.0035 

Weight 

0.06 

0.07 

0.11 

0.32 

0.06 

0.39 

- 

as(183 GeV) 

0.1111 

0.1076 

0.1117 

0.1032 

0.1075 

0.1084 

0.1079 

Statistical error 

±0.0036 

±0.0032 

±0.0032 

±0.0027 

±0.0035 

±0.0027 

±0.0027 

Experimental syst. 

±0.0027 

±0.0041 

±0.0033 

±0.0029 

±0.0028 

±0.0034 

±0.0031 

HERWIG hadr. corr. 

-0.0002 

±0.0020* 

-0.0009 

-0.0002 

-0.0005 

-0.0005* 

-0.0001 

ARIADNE hadr. corr. 

±0.0020* 

±0.0014 

±0.0014* 

±0.0006* 

±0.0020* 

-0.0003 

±0.0010* 

Hadronization error 

±0.0020 

±0.0020 

±0.0014 

±0.0006 

±0.0020 

±0.0005 

±0.0010 

Theory error 

±0.0041 

±0.0033 

±0.0049 

±0.0041 

±0.0042 

±0.0024 

±0.0036 

Weight 

0.15 

0.14 

0.10 

0.20 

0.12 

0.29 

- 



T 

Mu 


Bw 

c 

V23 

Weighted 

mean 

as(189 GeV) 

0.1146 

0.1070 

0.1121 

0.1026 

0.1079 

0.1066 

0.1075 

Statistical error 

±0.0020 

±0.0018 

±0.0019 

±0.0015 

±0.0019 

±0.0015 

±0.0016 

Experimental syst. 

±0.0019 

±0.0022 

±0.0014 

±0.0011 

±0.0015 

±0.0009 

±0.0012 

HERWIG hadr. corr. 

-0.0001 

±0.0022* 

-0.0009 

-0.0001 

-0.0005 

-0.0004* 

0.0000 

ARIADNE hadr. corr. 

±0.0019* 

±0.0014 

±0.0014* 

±0.0006* 

±0.0019* 

-0.0003 

±0.0007* 

Hadronization error 

±0.0019 

±0.0022 

±0.0014 

±0.0006 

±0.0019 

±0.0004 

±0.0007 

Theory error 

±0.0041 

±0.0032 

±0.0049 

±0.0040 

±0.0042 

±0.0023 

±0.0033 

Weight 

0.10 

0.13 

0.08 

0.15 

0.11 

0.42 

- 

aa(200 GeV) 

0.1142 

0.1048 

0.1128 

0.1030 

0.1082 

0.1064 

0.1068 

Statistical error 

±0.0023 

±0.0021 

±0.0022 

±0.0018 

±0.0023 

±0.0017 

±0.0019 

Experimental syst. 

±0.0024 

±0.0021 

±0.0034 

±0.0013 

±0.0028 

±0.0012 

±0.0010 

HERWIG hadr. corr. 

0.0000 

±0.0021* 

-0.0008 

-0.0001 

-0.0004 

-0.0003* 

±0.0001 

ARIADNE hadr. corr. 

±0.0018* 

±0.0013 

±0.0013* 

±0.0006* 

±0.0018* 

-0.0003 

±0.0006* 

Hadronization error 

±0.0018 

±0.0021 

±0.0013 

±0.0006 

±0.0018 

±0.0003 

±0.0006 

Theory error 

±0.0040 

±0.0032 

±0.0048 

±0.0040 

±0.0041 

±0.0023 

±0.0032 

Weight 

0.09 

0.15 

0.06 

0.16 

0.08 

0.45 

- 

as(206 GeV) 

0.1098 

0.1067 

0.1125 

0.1036 

0.1070 

0.1086 

0.1078 

Statistical error 

±0.0023 

±0.0020 

±0.0021 

±0.0017 

±0.0022 

±0.0017 

±0.0017 

Experimental syst. 

±0.0008 

±0.0019 

±0.0017 

±0.0009 

±0.0019 

±0.0015 

±0.0013 

HERWIG hadr. corr. 

0.0000 

±0.0021* 

-0.0007 

-0.0001 

-0.0003 

-0.0002 

±0.0001 

ARIADNE hadr. corr. 

±0.0019* 

±0.0013 

±0.0012* 

±0.0006* 

±0.0018* 

-0.0002* 

±0.0007* 

Hadronization error 

±0.0019 

±0.0021 

±0.0012 

±0.0006 

±0.0018 

±0.0002 

±0.0007 

Theory error 

±0.0040 

±0.0031 

±0.0048 

±0.0039 

±0.0041 

±0.0023 

±0.0033 

Weight 

0.11 

0.13 

0.09 

0.18 

0.09 

0.41 

- 


Table 10: Measurements of cts using event shape distributions at various values or ranges of c.m. energy. The values labelled 
189 GeV correspond to the data samples at 189 and 192 GeV, the values labelled 200 GeV combine the data at 196, 200 and 
202 GeV, and those labelled 206 GeV include all data above 202 GeV. The hadronization error is taken to be the larger of the 
effects observed using HERWIG and ARIADNE; in each case this is denoted by an asterisk. The weights and weighted mean 
are described in the text. 







T 

only 

Mh 

only 

C 

only 

B'y 

only 

Sw 

only 

y23 

only 

All 

ois{Mz) 

0.1242 

0.1181 

0.1177 

0.1222 

0.1134 

0.1193 

0.1191 

Stat. error 

±0.0011 

±0.0009 

±0.0011 

±0.0011 

±0.0009 

±0.0007 

±0.0005 

Expt. error 
Hadr. error 
Theory 

Syst. error 

±0.0018 

±0.0027 

±0.0053 

±0.0061 

±0.0014 

±0.0020 

±0.0042 

±0.0048 

±0.0015 

±0.0026 

±0.0053 

±0.0061 

±0.0017 

±0.0020 

±0.0062 

±0.0068 

±0.0015 

±0.0007 

±0.0051 

±0.0053 

±0.0014 

±0.0011 

±0.0031 

±0.0035 

±0.0010 

±0.0011 

±0.0044 

±0.0046 

Total error 

±0.0062 

±0.0049 

±0.0062 

±0.0069 

±0.0054 

±0.0036 

±0.0047 


Table 11: Combined as{Mz) fit results based on distributions of different observables, 
averaged over all c.m. energies, using OPAL data. 
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(Cl) 

(Up 




as(Mz)) 

0.1267 

0.1242 

0.1172 

0.1214 

0.1223 


Statistical error 

±0.0003 

±0.0002 

±0.0002 

±0.0003 

±0.0006 


Experimental syst. 

±0.0010 

±0.0008 

±0.0006 

±0.0013 

±0.0027 


HERWIG hadr. corr. 

-0.0017 

-0.0022 

-0.0023 

-0.0022* 

±0.0005* 


ARIADNE hadr. corr. 

±0.0038* 

±0.0035* 

±0.0029* 

±0.0010 

±0.0001 


Hadronization error 

±0.0038 

±0.0035 

±0.0029 

±0.0022 

±0.0005 


variation 

+ 0.0072 

- 0.0058 

+ 0.0066 

- 0.0053 

+ 0.0050 

- 0.0038 

+ 0.0050 

- 0.0004 

+ 0.0052 

- 0.0039 


X^/dof. 

i.2li 

6.2/3 

11.6/3 

2.2/3 

1.0/3 



((l-Tf) 

(cb 

(Hi) 

{Bh) 


(M|) 

as(Mz) 

0.1427 

0.1412 

0.1344 

0.1216 

0.1235 

0.1226 

Statistical error 

±0.0007 

±0.0004 

±0.0004 

±0.0006 

±0.0014 

±0.0003 

Experimental syst. 

±0.0017 

±0.0015 

±0.0012 

±0.0020 

±0.0031 

±0.0018 

HERWIG hadr. corr. 

±0.0006 

-0.0003 

-0.0011 

-0.0012* 

±0.0015* 

±0.0025* 

ARIADNE hadr. corr. 

±0.0041* 

±0.0040* 

±0.0029* 

±0.0002 

±0.0005 

±0.0020 

Hadronization error 

±0.0041 

±0.0040 

±0.0029 

±0.0012 

±0.0015 

±0.0025 

variation 

+ 0.0120 

- 0.0098 

+ 0.0118 

- 0.0096 

+ 0.0115 

- 0.0095 

+ 0.0050 

- 0.0036 

+ 0.0057 

- 0.0044 

+ 0.0045 

- 0.0030 

X^ldof. 

1.8/3 

2.5/3 

3.1/3 

1.3/3 

0.7/3 

0.8/3 


((l-Tp) 

(C3) 

(-®t) 

{Bi) 

{yg^) 

(Mg) 

aaiMz) 

0.1501 

0.1494 

0.1439 

0.1228 

0.1227 

0.1266 

Statistical error 

±0.0015 

±0.0006 

±0.0009 

±0.0014 

±0.0029 

±0.0006 

Experimental syst. 

±0.0020 

±0.0018 

±0.0017 

±0.0032 

±0.0040 

±0.0026 

HERWIG hadr. corr. 

±0.0018 

±0.0009 

-0.0002 

-0.0010* 

±0.0016* 

±0.0035* 

ARIADNE hadr. corr. 

±0.0044* 

±0.0043* 

±0.0028* 

-0.0002 

±0.0006 

±0.0015 

Hadronization error 

±0.0044 

±0.0043 

±0.0028 

±0.0010 

±0.0016 

±0.0035 

x^ variation 

+ 0.0144 

- 0.0117 

+ 0.0145 

- 0.0118 

+ 0.0149 

- 0.0121 

+ 0.0052 

- 0.0039 

+ 0,0056 

- 0.0042 

+ 0,0061 

- 0.0047 

X^/dof. 

1.0/3 

1.7/3 

1.5/3 

1.2/3 

0.4/3 

1.0/3 


((l-Tp) 

(cy 

(H|) 

{Bi) 

{yg^) 

(Mg) 

as(Mz) 

0.1557 

0.1553 

0.1495 

0.1222 

0.1214 

0.1276 

Statistical error 

±0.0028 

±0.0007 

±0.0017 

±0.0030 

±0.0057 

±0.0009 

Experimental syst. 

±0.0021 

±0.0019 

±0.0020 

±0.0051 

±0.0085 

±0.0034 

HERWIG hadr. corr. 

±0.0026 

±0.0016 

±0.0004 

-0.0011* 

±0.0015* 

±0.0040* 

ARIADNE hadr. corr. 

±0.0050* 

±0.0048* 

±0.0027* 

-0.0005 

±0.0006 

±0.0010 

Hadronization error 

±0.0050 

±0.0048 

±0.0027 

±0.0011 

±0.0015 

±0.0040 

x^ variation 

+ 0.0164 

- 0.0132 

+ 0.0168 

- 0.0134 

+ 0.0173 

- 0.0139 

+ 0.0049 

- 0.0035 

+ 0.0053 

- 0.0040 

+ 0.0063 

- 0.0049 

X^ldof. 

0.6/3 

1.3/3 

0.9/3 

1.0/3 

0.3/3 

1.2/3 


{(l-Tf) 

(C5) 

(H|) 

(sk) 

{yg^) 

(Mg) 

as(Mz) 

0.1615 

0.1606 

0.1541 

0.1217 

0.1196 

0.1269 

Statistical error 

±0.0053 

±0.0009 

±0.0033 

±0.0063 

±0.0111 

±0.0013 

Experimental syst. 

±0.0027 

±0.0019 

±0.0026 

±0.0086 

±0.0207 

±0.0040 

HERWIG hadr. corr. 

±0.0032 

±0.0022 

±0.0007 

-0.0011* 

±0.0013* 

±0.0042* 

ARIADNE hadr. corr. 

±0.0058* 

±0.0052* 

±0.0026* 

-0.0008 

±0.0006 

±0.0007 

Hadronization error 

±0.0058 

±0.0052 

±0.0026 

±0.0011 

±0.0013 

±0.0042 

x^ variation 

+ 0.0185 

- 0.0147 

+ 0.0189 

- 0.0150 

+ 0.0195 

- 0.0154 

+ 0.0044 

- 0.0030 

+ 0.0049 

- 0.0037 

+ 0.0062 

- 0.0047 

X^/dof. 

0.4/3 

1.0/3 

0.5/3 

0.7/3 

0.2/3 

1.3/3 


Table 12: Measurements of as{Mz) from event shape moments over the full range of 
c.m. energy, 91-209 GeV. The hadronization error is taken to be the larger of the effects 
observed using HERWIG and ARIADNE; in each case this is denoted by an asterisk. 
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(1-T) 


Figure 1; Distributions of thrust, (1 — T), at four c.m. energy points — 91 GeV, 133 GeV, 
161-183 GeV (labelled 177 GeV) and 189-209 GeV (labelled 197 GeV). The latter three 
have been multiplied by factors 3, 9 and 27 respectively for the sake of clarity. The inner 
error bars show the statistical errors, while the total errors are indicated by the outer error 
bars. The predictions of the PYTHIA, HERWIG and ARIADNE Monte Garlo models as 
described in the text are indicated by curves. The lower panels of the hgure show the 
differences between data and Monte Garlo, divided by the total errors, at 91 and 197 GeV. 
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0 0.1 0.2 0.3 0.4 0.5 0.6 

Mh 


Figure 2: Distributions of heavy jet mass, Mh, at four c.m. energy points — 91 GeV, 
133 GeV, 161-183 GeV (labelled 177 GeV) and 189-209 GeV (labelled 197 GeV). The 
latter three have been multiplied by factors 3, 9 and 27 respectively for the sake of clarity. 
The inner error bars show the statistical errors, while the total errors are indicated by 
the outer error bars. The predictions of the PYTHIA, HERWIG and ARIADNE Monte 
Garlo models as described in the text are indicated by curves. The lower panels of the 
hgure show the differences between data and Monte Garlo, divided by the total errors, at 
91 and 197 GeV. 
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Figure 3: Distributions of the C-parameter at four c.m. energy points — 91 GeV, 133 GeV, 
161-183 GeV (labelled 177 GeV) and 189-209 GeV (labelled 197 GeV). The latter three 
have been multiplied by factors 3, 9 and 27 respectively for the sake of clarity. The inner 
error bars show the statistical errors, while the total errors are indicated by the outer error 
bars. The predictions of the PYTHIA, HERWIG and ARIADNE Monte Garlo models as 
described in the text are indicated by curves. The lower panels of the hgure show the 
differences between data and Monte Garlo, divided by the total errors, at 91 and 197 GeV. 
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Figure 4: Distributions of the total jet broadening, i?T, at four c.m. energy points — 
91 GeV, 133 GeV, 161-183 GeV (labelled 177 GeV) and 189-209 GeV (labelled 197 GeV). 
The latter three have been multiplied by factors 3, 9 and 27 respectively for the sake of 
clarity. The inner error bars show the statistical errors, while the total errors are indicated 
by the outer error bars. The predictions of the PYTHIA, HERWIG and ARIADNE Monte 
Garlo models as described in the text are indicated by curves. The lower panels of the 
hgure show the differences between data and Monte Garlo, divided by the total errors, at 
91 and 197 GeV. 
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Figure 5: Distributions of the wide jet broadening, i?w, at four c.m. energy points — 
91 GeV, 133 GeV, 161-183 GeV (labelled 177 GeV) and 189-209 GeV (labelled 197 GeV). 
The latter three have been multiplied by factors 3, 9 and 27 respectively for the sake of 
clarity. The inner error bars show the statistical errors, while the total errors are indicated 
by the outer error bars. The predictions of the PYTHIA, HERWIG and ARIADNE Monte 
Garlo models as described in the text are indicated by curves. The lower panels of the 
hgure show the differences between data and Monte Garlo, divided by the total errors, at 
91 and 197 GeV. 
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Figure 6: Distributions of the two- to three-jet transition point, at four c.m. energy 
points — 91 GeV, 133 GeV, 161-183 GeV (labelled 177 GeV) and 189-209 GeV (labelled 
197 GeV). The latter three have been multiplied by factors 3, 9 and 27 respectively for 
the sake of clarity. The inner error bars show the statistical errors, while the total errors 
are indicated by the outer error bars. The predictions of the PYTHIA, HERWIG and 
ARIADNE Monte Garlo models as described in the text are indicated by curves. The 
lower panels of the hgure show the differences between data and Monte Garlo, divided by 
the total errors, at 91 and 197 GeV. 
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Figure 7: Distributions of thrust major, Tmaj., at four c.m. energy points — 91 GeV, 
133 GeV, 161-183 GeV (labelled 177 GeV) and 189-209 GeV (labelled 197 GeV). The 
latter three have been multiplied by factors 3, 9 and 27 respectively for the sake of clarity. 
The inner error bars show the statistical errors, while the total errors are indicated by 
the outer error bars. The predictions of the PYTHIA, HERWIG and ARIADNE Monte 
Garlo models as described in the text are indicated by curves. The lower panels of the 
hgure show the differences between data and Monte Garlo, divided by the total errors, at 
91 and 197 GeV. 
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Figure 8: Distributions of thrust minor, at four c.m. energy points — 91 GeV, 

133 GeV, 161-183 GeV (labelled 177 GeV) and 189-209 GeV (labelled 197 GeV). The 
latter three have been multiplied by factors 3, 9 and 27 respectively for the sake of clarity. 
The inner error bars show the statistical errors, while the total errors are indicated by 
the outer error bars. The predictions of the PYTHIA, HERWIG and ARIADNE Monte 
Garlo models as described in the text are indicated by curves. The lower panels of the 
hgure show the differences between data and Monte Garlo, divided by the total errors, at 
91 and 197 GeV. 
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Figure 9: Distributions of aplanarity, A, at four c.m. energy points — 91 GeV, 133 GeV, 
161-183 GeV (labelled 177 GeV) and 189-209 GeV (labelled 197 GeV). The latter three 
have been multiplied by factors 3, 9 and 27 respectively for the sake of clarity. The inner 
error bars show the statistical errors, while the total errors are indicated by the outer error 
bars. The predictions of the PYTHIA, HERWIG and ARIADNE Monte Garlo models as 
described in the text are indicated by curves. The lower panels of the hgure show the 
differences between data and Monte Garlo, divided by the total errors, at 91 and 197 GeV. 
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Figure 10: Distributions of the D-parameter at four c.m. energy points — 91 GeV, 
133 GeV, 161-183 GeV (labelled 177 GeV) and 189-209 GeV (labelled 197 GeV). The 
latter three have been multiplied by factors 3, 9 and 27 respectively for the sake of clarity. 
The inner error bars show the statistical errors, while the total errors are indicated by 
the outer error bars. The predictions of the PYTHIA, HERWIG and ARIADNE Monte 
Garlo models as described in the text are indicated by curves. The lower panels of the 
hgure show the differences between data and Monte Garlo, divided by the total errors, at 
91 and 197 GeV. 
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Figure 11; Distributions of sphericity, S, at four c.m. energy points — 91 GeV, 133 GeV, 
161-183 GeV (labelled 177 GeV) and 189-209 GeV (labelled 197 GeV). The latter three 
have been multiplied by factors 3, 9 and 27 respectively for the sake of clarity. The inner 
error bars show the statistical errors, while the total errors are indicated by the outer error 
bars. The predictions of the PYTHIA, HERWIG and ARIADNE Monte Garlo models as 
described in the text are indicated by curves. The lower panels of the hgure show the 
differences between data and Monte Garlo, divided by the total errors, at 91 and 197 GeV. 
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Figure 12: Distributions of oblateness, O, at four c.m. energy points — 91 GeV, 133 GeV, 
161-183 GeV (labelled 177 GeV) and 189-209 GeV (labelled 197 GeV). The latter three 
have been multiplied by factors 3, 9 and 27 respectively for the sake of clarity. The inner 
error bars show the statistical errors, while the total errors are indicated by the outer error 
bars. The predictions of the PYTHIA, HERWIG and ARIADNE Monte Garlo models as 
described in the text are indicated by curves. The lower panels of the hgure show the 
differences between data and Monte Garlo, divided by the total errors, at 91 and 197 GeV. 
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Figure 13: Distributions of the light jet mass, Ml, at four c.m. energy points — 91 GeV, 
133 GeV, 161-183 GeV (labelled 177 GeV) and 189-209 GeV (labelled 197 GeV). The 
latter three have been multiplied by factors 3, 9 and 27 respectively for the sake of clarity. 
The inner error bars show the statistical errors, while the total errors are indicated by 
the outer error bars. The predictions of the PYTHIA, HERWIG and ARIADNE Monte 
Garlo models as described in the text are indicated by curves. The fluctuations seen in 
the curves at low Ml are real artefacts of the models, while those at high Ml are merely 
caused by statistical fluctuations. The lower panels of the figure show the differences 
between data and Monte Garlo, divided by the total errors, at 91 and 197 GeV. 
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Figure 14: Distributions of narrow jet broadening, Bfq, at four c.m. energy points — 
91 GeV, 133 GeV, 161-183 GeV (labelled 177 GeV) and 189-209 GeV (labelled 197 GeV). 
The latter three have been multiplied by factors 3, 9 and 27 respectively for the sake of 
clarity. The inner error bars show the statistical errors, while the total errors are indicated 
by the outer error bars. The predictions of the PYTHIA, HERWIG and ARIADNE Monte 
Garlo models as described in the text are indicated by curves. The lower panels of the 
hgure show the differences between data and Monte Garlo, divided by the total errors, at 
91 and 197 GeV. 
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Figure 15: Moments of thrust, ((1 — T)""), heavy jet mass, (Mg), and C-parameter, (C”), 
for n = 1,..., 5 at four c.m. energy points — 91 GeV, 133 GeV, 161-183 GeV (labelled 
177 GeV) and 189-209 GeV (labelled 197 GeV). The inner error bars show the statistical 
errors, while the total errors are indicated by the outer error bars. The predictions of 
the PYTHIA, HERWIG and ARIADNE Monte Garlo models as described in the text are 
indicated by curves. The lower panel of the hgure shows the differences between data and 
Monte Garlo, divided by the total errors, as a function of n for the 91 GeV and 197 GeV 
data. 
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Figure 16: Moments of total jet broadening, (-Bx); wide jet broadening, (B^), and 
Durham jet resolution parameter, {(|/^)"'), for n = 1,..., 5 at four c.m. energy points — 
91 GeV, 133 GeV, 161-183 GeV (labelled 177 GeV) and 189-209 GeV (labelled 197 GeV). 
The inner error bars show the statistical errors, while the total errors are indicated by 
the outer error bars. The predictions of the PYTHIA, HERWIG and ARIADNE Monte 
Garlo models as described in the text are indicated by curves. The lower panel of the 
hgure shows the differences between data and Monte Garlo, divided by the total errors, 
as a function of n for the 91 GeV and 197 GeV data. 
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Figure 17: Moments of thrust major, ), thrust minor, ), and sphericity, (S”), 
for n = 1,..., 5 at four c.m. energy points — 91 GeV, 133 GeV, 161-183 GeV (labelled 
177 GeV) and 189-209 GeV (labelled 197 GeV). The inner error bars show the statistical 
errors, while the total errors are indicated by the outer error bars. The predictions of 
the PYTHIA, HERWIG and ARIADNE Monte Garlo models as described in the text are 
indicated by curves. The lower panel of the hgure shows the differences between data and 
Monte Garlo, divided by the total errors, as a function of n for the 91 GeV and 197 GeV 
data. 
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Figure 18: Moments of oblateness, (O'^), light jet mass, (Ml), and narrow jet broadening, 
(S^), for n = at four c.m. energy points — 91 GeV, 133 GeV, 161--183 GeV 

(labelled 177 GeV) and 189-209 GeV (labelled 197 GeV). The inner error bars show 
the statistical errors, while the total errors are indicated by the outer error bars. The 
predictions of the PYTHIA, HERWIG and ARIADNE Monte Garlo models as described 
in the text are indicated by curves. The lower panel of the hgure shows the differences 
between data and Monte Garlo, divided by the total errors, as a function of n for the 
91 GeV and 197 GeV data. 
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Figure 19: Fits of ttg to event shape distributions for 1 — T, Mh and C at 91 GeV 
and 197 GeV (i.e. combined 189-209 GeV). Each data point shows the measured bin 
contents divided by the integral of the predicted distribntion across the bin; the inner error 
bars indicate statistical nncertainties, and the onter bars show the combined statistical 
and experimental contribntions. The dashed curves represent fractional variations in the 
predicted distribntions, corresponding to our perturbative theory uncertainties in ag. The 
slightly wider shaded bands indicate the combined theory and hadronization nncertainties. 
The ranges used for htting each distribution are shown by horizontal arrows. 
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Figure 20: Fits of ttg to event shape distributions for i?T, -Bw and at 91 GeV and 
197 GeV (i.e. combined 189-209 GeV). Each data point shows the measured bin contents 
divided by the integral of the predicted distribution across the bin; the inner error bars 
indicate statistical uncertainties, and the outer bars show the combined statistical and 
experimental contributions. The dashed curves represent fractional variations in the pre¬ 
dicted distributions, corresponding to our perturbative theory uncertainties in cts. The 
slightly wider shaded bands indicate the combined theory and hadronization uncertainties. 
The ranges used for htting each distribution are shown by horizontal arrows. 
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Figure 21; Measurements of cts using fits to distributions of six event shape observables. 
The inner error bars represent statistical uncertainties and the outer error bars the total 
uncertainties. The grey symbols indicate, without errors, previously published OPAL 
measurements, which are superseded by our new results. 
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Figure 22: OPAL combinations of as{Mz) values inferred from distributions of individual 
event shape observables. The inner error bars are statistical, while the outer bars represent 
total uncertainties. The grey band corresponds to the total uncertainty of the combined 
as{Mz) value, and the dashed lines indicate its statistical uncertainty. 
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Figure 23: OPAL as{Mz) combinations based on distributions of individual event shape 
observables at different energies. The inner error bars are statistical, while the outer bars 
represent total uncertainties. The grey band corresponds to the total uncertainty of the 
combined a^^Mi) value, and the dashed lines indicate its statistical uncertainty. 
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Figure 24: A global QCD running fit to the OPAL as measurements based on event 
shape distributions. Each point represents a fit to the six measurements at an individual 
centre-of-mass energy, while the curve represents a global £t to all measurements. The 
form of the curve is determined by the 0{a^) Renormalization Group Equation of QCD, 
with as{Mz) as a free parameter. The grey band corresponds to the total uncertainty of 
the fitted as(Mz) value, and the dashed curves indicate the statistical uncertainty. 
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Figure 25: Measurements of Og using fits to moments of six event shape observables. The 
inner error bars represent statistical uncertainties, and the outer error bars show the total 
errors. The dashed line indicates the weighted average described in the text; only the 
measurements indicated by solid symbols were used for this purpose. 
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